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DEVELOPING THE FUNDAMENTAL THEOREM OF CALCULUS 



1^ INTRODUCTION 



One of the jpost useful tools of matliematics is the 
Fundamental Theorem of Calculus^ Although beginning 
calculus students learn to use the Fundamental Theorem of 
Calculus to' find the value of a definite integral, they 
are often ieft with only a manipulative tool and not a 
thorough understanding of the tool itself, its proof, and 
Its numerous applications. The purpose of this module is 
to develop, by means of sOme applications, some of the 
basic concepts of the Fundamental Theorem of Calculus. 



2.# THREE SIMILAR PROBLEMS 

ConsideV'the three fol*owi«g. -problems , followed by 
their solutions : • 

^ * f ' *^ 

EXAMPLE 1. Find the area of the rectangle that li 15 feet 
by 10 feet. 

EXAMPLE 2. Find the ^distance traveled by a car if its 

velocity is 10 feet per second and the length 
of time it travels i^s IS seconds. 

EXAMPLE 3. Find t^e work required to lift a 10-pound bag 
of salt a distance of 15 feet. 

SOLUTIONS. All th(-ee .problems have common numerical quan- 
tities, namely the 10 units and the 15 units. 
Figure 1 and Table 1 show^the relationship 
between these quantities and the solution of 
each example. 

It should be noted in each pf thes^ examples that the 
vertical quantity is a constant function of the horizontal 
quantity, ' i. e. , in EXAMPLE 1, the width w is%lways 10 feet 



for 
for 
0 < 



0 < i < 15, in EXAMPLE 2; the' Velocity .y = 10 ft/sec 
0 < t < 15, and in EX.\MPLE 5, the force 'f = 10 lb *for 
d < -15. ' ^ 




* 15 

f^'9ure r; Solution to EXAMPLES i; 2, and 3. 



TABir I P 
Solution to EXAMPLES 1, 2, and 3.> 



EXAM- 
PLE ' 


Hoi^i zonta 1 
Axis 


Vertical 
Axis ' 


Formula 




. .Sol ut ion 




I 


- 15 ft 


w 


= 10 ft 


A = i xw 


Area 


= 15x10 = 


150 ft^- ^ 


2 


t = 1$ sec 


V 


= JO "ft/sec ' 


0 = V X t 


Dist 


.= 15x10 = 


150 ft 


3 


d = 15 Vt 


f 


= 10 lb 


W = f xd 


Work 


= l'5x fo = 


150 ft|-lb 



r 



Exerc ises 

1. A car, nx^ying with a velocity of 10 ft/sec begins to^ reduce its 
speed uniformly* and comes to rest in 15 seconds. Explain why the 
product of the velocity and, time does not produce'the distance it 
CT^/e Is. 

f 

2. A 1,0-pound bag of salt is lifted up steadily to a height of J5 feet. 
'The bag has a hole In it and salt leaks- ou^ at a uniform rate, so 

* ithat the bag is empty when it r<eaches the destined height. Explain 
why the product of tfie force and the distance does not produce the 
v/ork necessary to lift the bag. 



niRHH MORE SliMILAR P^LEMS 



5,1 St-a tement of the Problems 



Consider t}>e •-th re e* £ol lowing problemsi ' 

FrTTd^t'he^ ,area of a' r»t§ht triangle whose base 
is 15 feet and whose altitude is 10 feet. 



EXAMPLH 4 



EXAMPLE 6'! 
9 



EXAMPLE S. Find >the distance traveled by a car whose 

velocity at time t xs v(t) = ((-2/5)t + 10]' 
ft/sec^ when It travels fVom t = 0 to t = 15* 
'seconds. (See Exercise 1.) ^ * 

Find the work done in lifting a^bag of sal.t>a 
distance of 15 feet ab.ove the ground, .assuming*^ 
that the bag has a hole in it so fhal^.at height 
' ^d above the ground its weight (m^^nitude of 

downward ^orce) i^ g(d) = ((-2/3)d + 10] lb, ^ 
(See sExercise' 2.) . , 

5.2 *^olution of One of the Probl.ems ^ '> 

Ail three^examples (EXAMPLES 4, 5, and 6) have been 
purposely set up to jjialce -their solutions - sinhiiar . .(Thi^ 
mjght no_t be noticeab.4e at a fiTst glance!) For EXAMPLE 4, 
consider Figure 2 to you in its solution^ 




15 % ' 

Figure 2. Graphical aid /or solution to EXAHPLC k. 



ERIC 



Since the area of a triangle =^ [ 1/ 2^) [ 1 ength o.f babe) (length 
.0/ altitude), then^ the area * ( 1/ J) i, 1 S Ji 1 Q) 75 square 
feet. It bhould_J)e remembered that 'this^ problem was solved 
rapidly by using an* appropriate formuja. , - 

Notice that in Figure *2, w 'is jiot a, constant functicfn 
of i. Since the line segment 1% Figure 2 Kas'^ slope ef 
^(10-0)/ (0- 15) = (-2/»5) and a w-intercep^ of 10,/ then 
W(i) = (-2/3)1 + 10^, where 0 < i < IS. For the functions 
m EXAMPLES 5 a'nd 6, their graphs ar^ given in Figures 5 
and 4 respecpvjely , where v(t) = (-2/5)t + 10 for 0 t 15 
and g(d) = (-2/5)d + 10 for 0 < d < 15. 




15 t . 15 d 

, Figure 3- Graph of function v. Figure*^. Graph of function g. 



The fog th*a.t may'-have hung ov.er EXAMPLES 4, 5, ani 6 
should' now be clearing and one should begin to see the , 
similar ity«. between these three examples'. *Tor all three 
examples, we have the same function and the same doMin, 
but each with different labels: 

/. EXAMPLE 4: w(Jl) =^(-2/3)£ + 10, where 0 < £ ^ 15. 
EXAMPLE 5: v(t) = (-2/3)t + 10, where 0 < t < 15. 
EXAMPLE 6: *g(d) = (-2/3)d + ;0, where 0 < d < 15r 

qUESTION: Since the solution to EXAMPLE 4 is 75 s/^uare 
feet and each example has the same' function (with different 
labels), c^n we assume that EXAMPLES 5 and 6 will have the 
same num'erical solution but with different label's? 



To answer this question, we wiu'go back'; and solve 
EXAMPLE 4 by a different method, where wc do, not^ use the 
area of-a triangle formula. This will enable us to also 
solve the two other examples. 

5.5 Approximate Solutions of EXAMPLES 4 , . 5^, and 6 

0 -We will start by first approximating the area of the 
triangle in Figu-re 2 by a sum of areas of rectaA^les. 
First we divide (0,15) into 15 subintervals of equal length; 
each of thfe form Ii-l,i), where i»= 1, 2, 5, . . . , 15 as- 
shown in Figure 5. Over th^ ith subinterval, (i-l,i), we 

4 



-1 L. 



-1 L 



° 1 .2 ^ 5 6 7 8 9 10 II 12 13 1^1 15 

Figure 5. Subd i vi s ion' of [Q, 



A 



construct the ith rectangle whose "length" is the length 
of the Ith subinterval and whose "width" is determined by 
the right endpoint 1, as %hown m Figure 6, » The area of 



1 St 

\ rectangle 
,2".^ rectangle 




.th 

I reptangle 



Figure^.' Approximation of area of'^triangfJe. 



the Ith rectangle is the product of the width w(i) and the 
length of the ^th suljinterval , namely i (i - 1)-= 1, or 



(A) , w(i).J = [(.2/3)i^+ 10] (1) \ ' 

for i = 1, ^2, 3, . . . , 15. Since there are '15 such Rec- 
tangles, then the area'of all of these is - i 

• r5 15 ' ' ' 

(B) I wtn-l = I ((-^/3)i + 10] (1), as w(i) = (-2/3)i +'lO, 



i = l 




i = l 

15 15 « 

= I (-2/^3)1 > I (lO)-l , 
1 = 1 X -1 = 1 

' 15 ' * 15 . 
= .1-2/3) I 1 ^ (10) J 1, 
i = l • . i = l 

= (-2/3 

v( 1 5 ) ( 1 5 + 1 ) ^ 10(1.5), 



s I i = nli^'and I 1 = rf. 



j. = l 



1 = 1 



= 70. 

. \ ^ 

Hence, the area of the triangle is approximately 70 square 

feet. \ ' ' \ 

Let's'Vo bacl^ and Iodic Vt- (A) , (B) , and (C) " in t^rms 
*of EXAMPLE and Figures 3 and %, Changing the labels, we 
have the apntoximate- distance over thq ith subinterval 
represented-^y 

^(A/) v(i:)Vl - I(-2/3)i + ip](l), 

whete v(ij = ('fi-2/3)i + IQ] is the velocity over the i th 
subinterval and\ r represents 1 second ofx elapsed time. 
Notice that we ^re assuming that the velocity is a cons tan^t 
value, ovfer the ith subinterval.' In turn. 



(B») 



15' I . 
I v(i)-l J 
i = l ' * 



IS an approximation to the distance traveled over* the 15- 
seco-nd time perio^. This ^um worlcs out to be 



< 



(K*) I v(i)-l = 70 ft. 

i = l 



,. 2- 



Exerc i ses ' \ • ^ 

3- State what (A), (B) , and (C) represent with* the appropriate labels 
when the width function w is replaced, by weight function g of 
Example 6. . 

^ k. Approximate the area of the triangle ih EX>AMPi.E k by divitiirig the 
interval lOJSr^'into ^5 subintervals. Use che left endpoint 
to determine the height of each rectangle. 
5a. Interpret/^the results of Exercise above in termS of EXAMPLE 5. 
.b. interpret, the results of Exercise k ibove in terms of EXAMPLE 6 

'■ • - • : - 



\. ' o Since .70 (with the'^proper label) is an approximation 
for the solution of EXA^iPLES 4, 5, and 6, iets so^ve all 
of these problenvs by a more general nfethod. Let 

= (-2/3)x + K^>whcre 0 f ^ £ 15, "where function f 

, represents^any .one of. the three functions in EXAMPLES 4, 5. 

,and|6. Generai i z mg Figure 5^ we divide I0,'l5} into n 
equal length subrntervals (each^ofjength 15/n) and denote 
the ith subintervaO^ by l( ^nfT^/^; ( i) ( 1 S/q.) ] / for 
1 = 1 2 , 3 , . . / ^ n as show^i m Figure*?. 



15 



Figure ;. -.Subdivision of«[0. I5) into n equal length subinterval^ 



Consider then the product 
(D) - ft(i)(15/n)],(l5/n), 



rwhicVmay be interpretecJ in the fallowing three ways: 

1. For EXAMPLE 4, (D) is the area of the.ith rectangle 

Vith width f[U)115/n)) and length (15/n), where the 

width IS determined by the right* endpoint , (i)(15yn), 

of the i th 'subinterval . 

•J ' . 

2. For EXAMPLE 5, (D) is the approximate distance 

^ traveled over the i_th subinterval of time length (IS^n) 
and where the velocity, f ( ( i ) ( 1 5/n) ) , i» constant^ over 
the interval and it i^ determined by the Vight endpoint , 
- (i){15/n)', of the itji subinterval. - 

3. For EXAiMPLL 6 ,\ (;D)'is»'the approximi^te work done over 
the ith sub' interval of distance (15/n) and where the * 

^ force, f ( ( 1 ) ( 1 5/n) ) , is cpns tant over the interval a«d 

it IS determin^(^ bj* th*e right erfdpoint ,^ ( i )^(1 5/n) , of 
' the ith subljjterval ^ ^ ^ J 

Adding up all n of these^ product ,* we then have,^ ■ ' 

' . ^ \ ^ * 

(E) ' I f'UiHlS/n))(15/n) I 2/5) ( i ) (1 5/n) + 10) (15/n) 
i=i . i=l 



\ [(-loI/nO ^+ 10] ClS/n), 
1 = 1 



n 



((-15^?^2) + l50/n] , 
i=l ' ' ^ 



t 



Vj50/n^) v\ (i)]. 
if^l 



^ (150/n) ( I (f 
\ (^l'50/n^)(i^^i-^^).^^. (150/Ti)(n) 



or 



n 



(F) I f 1(1) (i5/n})(15/n) = (-75^(1 + l/»n) + 150. ^. 

1=1 * ^ ^ ^ ^ , 



Notice that the expression in (F) is a function only 
of a, the number of submtervals of [0,15). By letting n^^ 
take on specific values in (F), we obtain various approxi- 
.mations to the solutions of EXAMPLES 4, 5, and 6, For the 
case thatvn * 15, we obtain the numerical valae of 70, the* 
result we had previously seen. . ' » 

• ' Returning, to (F), we will ^et n->?o,* so that '* S 
(G) liin I f[(i){lS/n})C13/K) ^ ■ lim [ ( - 7 5) ( 1 + l/n) + 150), 

(-75)(1) + 150, ^ 

f 

=^ 75. 

Notice tKat the 75 is the same numerical value that we had 
fireviously obtained (see page 4) for the area of the tri- 
angle when we t^ad used the specific formula for determining 
the area of a triangle. Notice also, that we have now', 
answered ^he- question that was' posed on page 4. ' 



Exerci sesf 



6a. 



Find the value of the expression in (F) for the case that ' 
n = 25; n = 75; n = 300. 

Wh7 is each successive value of rr a better approximation than 
each previous value of n? 

The sum of products in.(F) was determined by using' the r?^ht 
endpoint , (i)(l5/n), of each subinterval. 

Determine a similar sum of products by using the left endpoint . 
(i " 0(I5/n), of each subinterval. Simplify the result as much 
as- poss i ble. 

In your results of 7a above » let n = 25; n = 75; n - 300.. 

Iti your r^ults of 7a above, let'n->«. Then compare the res Jl t ' 
wi th (G) . ■ 



8. As in Exercise 7 above, 'J 

•■^ _ . 4'-- 

a. . Determine a sum of products by using* the midpoint of each sub- 

interval. Then simplify it as much as possible. 

b. In' your results of*8a above, let n = 25; n = 75; n '= 300. 

\ 

c. In your results of 8a above, let n-x» and compare the resJIt with 
* (G). • ' ^ 

1^ ' ^ 



ND "tmSi 



4. >RIEMAi\N SUMS AND TMINDEFINITE INTHGR^M 

V ' - * ' 

4.1 Def mjition' of Riemann>Sum ' ^ - 

We have been looking at some special case's of what is 
called a Riemann Sum and the definite integral . Let us 
now look at these more generaf concepts. 

Let f be a function defined on a closed interval (a,b) 

Let a = and b = and select (n-1) points x^, x^, x^, 

, . . , between x^ and ^j^o that a = Xq < x^ < X2 < 

. . • < X , <*x = b. These (n+1) points are said to 
' n"i^n ^- > 

partition the interval [a,b] into n submtervals , where 
the itji sUbmterval is denoted by [x^_j,x^), ^r 1 = 1, 2, 
3, . . . , n. These n subinterv^ils are said to form a 
Partition P .of [a,b] . See Figure, 'S. '- 



H 1— H-l— — M M- 



^2^3 '^-1 ^i • V2 Vl \ 
Figure 8. A partition P of [^,b] , 



Denoting the length of the i_th subinterval by Ax. , we havfe 
(H) AX. = - x^.^ 

For each value of i, sele^ one point c from the ith^ sub- 
interval, so that ^^^-^ 1 1 now*form the product 



10 



(I) , f(c J-ix^ , 

and thea the .^um of 'these n products, 



1 = 1 ' 



This sum of productS'(J) is called a Riemann Sum of t.he 
function f over ta,b] for t^e partition P a,nd the 
"choice of c^. It^ should be noted that a RiemAnn Sum 
depends upon 

1 . the function f ^ 

^2. the closed interval U,b] over which f is defined, 
5. the partition P of ta,b], where each submterval 

need not be of the same length, and 
4. the point c. selected from each' subinterval * 

Returning to the Riemann Sum in (E) , we see that 

1. the function f is given by f(x) = (-2/3)x + 10, 

2. [0,15], is the closed interval over which f is 
defined, 

^ 3. the partition P "of [0,1^5] consists of n sub- 
intervals q{ equal length, and ' . ' 
4. the point c. = (i)y5/n) is the right endpoint of 
each 3ubinterval. 



Exerci ses 

9. Set up and simplify a Riemann Sum for the function f(x) » 2x + ^ 
over [1,8], using a partition with subintervals Qf equal length 
and selecting the left endpoint of each subinterval for c". . 

\ln Exerci/Se 9 abov^, let a take on various positi>ve integer 

values. ^ 

\ 1 

Iruerpret the numerical- results of Exercises 9 and 1*0 above in 
teUs of EXAMPLES ^, 5, and 6. 



11 
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4.2 Definition of the. Def'inite Integral 

In a partition F of [a*,b), the length of the longest' 
•subinterval is called the norm of the partition p, and is 
denoted by ||P|! . F^r the Riemann Sum in (E) , since each 
subint'erval 15 of the same length, then |tP|| = (I'S/n).' 
We iiow come to one of the main topics of this module, the 
Riemano-r^-Dr -definite integral. 



Def ini tion : , (I) Let f be a function defined over [a,b). 
(2) Let P be a partition of [a,b) having n subinterval s , " 
where Ax. is the length of the it^ subinterval and c. is^. 
any point in the i_th^ sub i n terval . (3) If there exists a 
number L such that 



lim I f(c,)'Ax. L, 
iPll-O i = l ' ' ' 



then L is called the Riemann Integral or the definite 



integral of f over [a,b) and is denoted by f(x)dx; i.e., 



[ f (x)dx 



I im y f(c,) -Ax. 



P kO i = l 



Retorning to Riemann S\im (E) or (F) and the partition 
P associated wi^h this Riemann Sum, we have 1|P|| = (15/n). 
Since ||P||-^0 as n-^^, we then have (G) which now can be 
written as ^ 
fl5 



(K) 75 



0 



[(-2/3)x + lOjdx = lim I f[(i) (15/n)) (15/n). 
n-*'^ i-i 



Exerci ses 



12. In the results of Exercise 9 above, let n -x» to find the value of 
|^(2x + ^)dx. 



13a. Write a Riemann Sum to approximate the vajue (3x + 5)dx 

h 
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b. Find the vali;e of [^(3x + 5)dx by using the results of 13a above. 

Interpret thg-resOltS of 13b above, in terms o^ EXAMPLES A, 5, 
and 6. ' ' 

\ka. Find the value of - (5x + I8)dx. 

2 

b» . Find the value of ^(x + Jkx + 5)dx'. 



S> THE TUNDAMENTAL THEOREM OF XALCULUS 

QUESTION : Does one always have to evaluate a definite 
integral by calculating a limit of Riema'nn Sums? The 
answer to this question is. usually NO. For functions f 
th^t are, contiiiious over Ea,b] , the value of | f(x)dx can 
often, but not /always , be determined by the Fundamental 
Theorem of^Cal/ruIus . We will pre'sent this theorem in Sec-^ 
tion 5.3, but before we look at it we will look at some of 
its underlyii/g principles in EXAMPLES 7 and 8. 

S.l A Definite Integral with a Variable Endpoint 

In , we saw that [C-2*/3)x + lO]dx = 75. We now 
consider^ more general form of this integral, namely 

rx ^ - * ' - . * 

(L) / E(-2/3)t + 10]dt, 

•0 ^ ' . 

whex^ we will assume that 0 £ x' £ 15. Note that in (L), 
we/have changed f(x) = (-2/3)x + 10 to -f(t) = ('-2/3)t + 10; 
fad we consi4ered |q [(-2/3)x + 10]dx, then we would have 
'used X for two different purposes, namely, (1) to denote 
the right endpoint of EO,x], and (2) to denote the inde- 
pendent variable of the function f. 



Exercises ♦ v 

~ rx 



15a. If X = 0 in (L) , whatsis the value of 



[(-2/3)t + !0]dt7. 



0 

X 



i 



b. If X = 15 in (L), what is the value <Jf [ [(--2/3) t + IO]dt? 



.7 
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c. If X » 10 iti U) , what is the value of 
flO ^ 



'16. 



Interpret 



4(-2/3),t + I0]dt? ' , 



[(-2/3)t + lOjdt in terms of EXAMPLES 5, and 6. 



) 

J 



By letting \\ 0 1 ^ 1 15, take on various values in 
(L) , we obtain^one and only one defin^ite integral, so that^ 
we can consider the expression in (L) to be a function of 
X, namely 

(M) ' F(x) = KV2/5}t + lOJdt , 
where 0 _< x £ 15 and 0 _< t £ x . 

5.2 Antiderivat;ves ^and the Definite Integral ^ 

EXA>!PLE 7. Find a^nather representation of the function F 
that IS represented in (L)»and (^f) . • 

Solution : Since F(x)^is a definite integral with lovver 
limit 0 and upper limit x, we will find its value as- a 
limit of Riemann Sums. By dividing (0,x] into n equal 
length subintervals , then Ax. = x/n and the ith subinterval 
isM(i - l)(x/n),(i)(x/n)], for i = 1-, 2 , . , . , n , 
Selecting the right endpoint of each 'subinterval as C., 



then c. 
1 



(t)(x/n). See Figure 9, 




Figure 9. A partition of (0, x] . 



is , 
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For this partition P of,[0,x) and for f(t) = [(-2/5)t + 10), 
t-he Riemann Sum is 

n n * ' 

I f(c-)'Ax = I [(-2/5)(cJ + 10)(AxJ* 
1=1.- ' ^ i-1 1 ^ _ ^ 

n 

' ^ I K-Z/SJCi'x/n) + 10))(x/n), " 
1 = 1 

I [(-2x^73n-J(i) + (lOx/nHD) , 
i = l 



(-.2x'/3n') 



n 

i = l 





n 


+ (lOx/n) 


I (13 




1 = 1 



(NJ 



(-2x2/3n^)-[lllH^] . (10x/n)(n). 
2, 



= (-xV3)(l + 1/n) + lOx, ^ 

where thfe value of the Riemann Sum in (N) is a function of 
^Q^^ ^ (the right endpoint of [0,x)) and n (the number of 
equal lengthed subintervals of (0,x]). Since ||P|| = (x/n)->0 
as n-^-w, then 

F(x) = l(-2/3)t + IQJdt, 

io 



— = lim I f (c- ) -Ax. , 
llPlhO 1^1 ^ ^' 

= lim((-x^/3) (1 + 1/n) + lOx), 

2 .1 * ' 

• •= (-xV3)(l) + lOx, 

= (-x-^/3) + lOx. . . 

Hence, another* representation of the function F is 

(P) F(x) = (-x2/3) + lOx/ ' 

Upon investigating some of the properties o£ this 
function F in (P) and the given function f; we see that 



15 



'i) r'ix\ = 3lU-x-/3J + lOx) = l-2/5)x + 10, . 
2) fit) = (-':/3)U+ 10, or fix) = (•2/5)x + 10. 



Mere we see an extre^nely important concept , namely that 

'> I 

.f(^) = '(-x*'/5) + lOx IS an antiderivative of ^| 

.f(x) = (-2^3)x + 10; i.e., ' ' 



F' Cx) = f (xj ,'''"when F(x) 



[(-2/3)t + lOjdt. 



^0 



Note* Re'nember this fact when we investigate the Funda- 
mental Theorem of Calculus in Section 5.3. 



Exe rcises 

17. ' As in EXAMPLE* 7, f^nd another representation of the function F 

in (m) by letting c. = (i - l)(x/n), the left endpoint of each 
subthterval. Compare this result with (P). 

18. As in EXAMPLE 7, find another representation of the function F 
in (M) by the following method: 

a. In Figure 9, let G(x) = the area of the trapezoid whose two 
, parallel sides are 10 and (-2/3)x + 10 units \n length and whose 
base is x units In length. 

I?. Use the area formula for a trapezoid to determine G(x). 

Compare G{x) with the results of (p) in EXAMPLE 7- " ■ \ 

19. As in EXAMPLE 7, give another justification of the function F 
in (p) on the bast's that velocity is the derivative of a posi- 
tion function. 

Before looking at the statement of the Fundamental 
Theorem of Cal'culus, we will consider another example. 

EXAMPLE 8. Find the distance traveled from t = 2 to t = 8 

* 

seconds by a car whose velocity is 
v(t) = (-2/3)t + 10 ft/sec. 
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Solut.ion ; As *a definite iYitegral, this distance can be 
represented as l(-2/5)t + lOJdt; (Take a moment to con- 
sider why this is so.) Since r(x).= ((-2/5)t + 10]dt is 
,th€x ^distance the car travels dyer (O^x], then 

r8 ' ' ' . ' ' 

^3 F(8) ^ j- l(-2/5)t + lOldt IS the distjance it 



I) 



travels over [0,8] , and 



^H2) = l(-2/3)t + lOjdt IS the distance it 



travels over IP,2] . 

Since we are concerned vvith the distance travelecf over 

12,8), we then seek F(8), - F(2). That is, ^ ' 

!^ ' y " ^ ' 

l(-2/5)t + lOldt = F(g) - F(2i. 
But m CP) we saw that r(x) = (-xV3) + lOx, so that 

j^[(-2/3)t + lOjdt = [(-8V5)+10(8)],;- ((-2^/3) + 10(2)] ; 



= 40 feet. . ' ' 

5>5 Statement o f the Fundamental Theorem of Calctrlus 

We are now ready to state the Fundawiental Theorem' of-^ 
Calculus. 



FUNDAMENfAl THEOREM OF CALCULUS 

If (1) function T is continuous over (a,b] , and • ' 
(2) F(x) is an;^ antiderivati ve of f(x) over [a,b] , 
then . . < '. * 

- ■ f f(x)dx = F(b) - F(a). 

a. 
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^Returning to EXAMPLE 8, we see ,t^af\C^hav^e act^lly * 
use*&*the Fundamental Theorem of Ca Ic^'tni^ s^irj^ -tTtc $<>Iuti'on 
of the* problem, -as* t . , 'L^'^'^V"*^^ ' 'V/^^^ 

1) v(t). = (-2/5)t + 10 (or, ^^f^^^^^^/^^^'i^ 

^cantmuous over (2,8) -^s f uncti^^*^^CX.f;5n^^on 



f) 13 "a linear • funct-ion. ' . T^S?^'!^ l^*" 

^ ,2) F(t) .= (-t^/s) V i(rt (or„F(x3 = CVxV%:^^,^ 

IS an antiderivative of v(t) = ( - 2/>)*^'«»+^0tt|;tjW:i^.- - 
f(x)''= r-2/3)x ^ 10) over (2,8], - af^^ , , ./^i," C 

Hence, the two parts of the hypothesis of the xtieb/^'JVVye* ' 



been satisfied. In conclusion then, ' 



(Q)^ J^'[(-2/5)t + lOjdt - Ft8) - F(2j = 

the same result that we hadf previously obtained, c> r 



Exerci ses 

20. In the second hypothesis of the Fundamental Theorem of Cajculus, 
* it is sta'ted that F(x) is any antiderivative of f\x) o\f^r\[^ ,by. 

2 • ' , ^ 

a, gvaVuate (Q) using (-t /3 + lOt + 5) as an ant icter i vat i ve.^ ? 

b* Evaluate (Q) using (-t^/3 + iOt - V3J as an antider^tfve. 

c. Evaluate (Q) using (-t^/3 + IOt + k) as an ant ideri vat i.ve , 
where k is any constant. 

21. Return to item,(K) an^ Exercise aod evaluate each of tK§se 

definite integrajs by using the Fundamental .Theorem of Calculus. 

I * 

22. Use* the Fundamental Theorem of Calculus to evaluate each K)f fhe 
following definite integrals: ^ 

r5 



(x^ + 2x + l)dx 



2 

[ (x^)dx + li (x)dx +[ (\)d> 



d. ^ 



(x + I)dx> 



■J (x + I)dx. 



23- Explain why the Fundamental Theorem of Calculus cannot be used 
to evaluate the folloviing integraJ: •» ♦ 

' f (l/x^)dx. 

WARNING: Do not^ attemj^t to use the Fundamental Theorem of 
Calculus unles^all conditions of the hypothesis have been 



satisfied. 



2k. Evaluate each of the following by using" the Fundamental Theorem 



of Calculus : 



f (-3x^ + 5x - 7)dx 
J^'^(l/x^+ x"^)dx. 



b. (x + l/x^)dx' 



d. 



I^Kx^ +J)/x^]dx. 



25. Find the area of the regioa bound by the x-axis and the 
parabpla y = 6 - - x . 

26. Let v(t) = t^ - 3t + 2 be the velocity function of a^clr when 
0 £ t <^*^/where the velocity is measured ir^ ft/sec and the 

* time t is measured in terms of seconds. Find *the distance it 
travels. WARNING: What does a negative velocity indicate? 

27. A bag of salt originally weighing 1^^ pounds^ is lifted upward. 
The salt leaks out uniformly at a rate so that half of the salt^ 

V* is lost when the bag has been 1 i f Jed 18 feet. Find the work 
done in lifting the-bag this distance* 
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CONCLUSION 



The FundamentJvl Theorem of Calculus "has ^ long and 
.fascinating' history behind it... Prior to its de^velopment;*' 
mathematic^ians worked for centuries with the derivative, 
the anj:i^rivative, and sums of products.! Isaac Barrow 
tl620<fi^77) , a teache.r of Isaac Newto-n (1642- 1 727), di's- 
covered and proved the Fundamental Theorem of Calculus, 
although his method and terms w5re quite different from 
those used in this module. Using Barrow's results, both 
Newton and Gottfried Leibnitz (1646-.1716J , working inde- 
pendently of each otheV, developed many of the concepts^ 
of calculus, although much the calculus that we know 
and use today is attributed to Georg B. Riemann (1826- 
1866J. ^ 

/ * . 

In this module, we have solved area, distance, a'nd 
\iOx\ problems by the * Fundamental Theorem of Calculus. It 
can also bq used to.solve problems dealing with the .volume 
of a solid of revolution»^arc length, moments, center of 
gravity, hydrostatic force, product cost, growth (or decay) 
of a substance, etc. Good luck on your usa^e of this 
remarkable tool I ' ^ 

7. ANSWERS TO EXgRClSES. ' 
\ 

1. The formula, distance - velocity x time, can be used c^Iy when 
the velocity is a constant* which is not the casfe herej 

2. The formula, work = force x distance^ can be used only when the 
force is a constant, which is not the ca^ here. 

3- (A") g'(i)-l = [(-2/3)i + I0](1) is the approximate work done 

over the ith subinterval, where g(i) is the constant force over 

itih subinteryal and 1 represent^ 1 foot, the 'length of the 
i th subi nterva"! . 

\ 

/B")* 15 ^ ' 

\ g(i)*I represents an approximation of the work done in y 
i=l ' - . " ' * . ^ 

•»^»ng the bag of salt a distance of 15 feet. 
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(C") 70 represents the numerical approximation of the work 
done in terms of foot-pounds. 

k. The length of each subinterval as (15/^5) = 1/3 units, so that 

the left endpoint is i\ - 1) (1/3) , for i = I , 2; 3 ^5- 

The area of the i_th rectangle is w((i - I ) (1/3) ) (I /3)'. The 
.area of the kS rectangles is . ^' 

^5 ^ 

* I w((i - l)(l/3)) (r/3) = 76 2/3 sq. ft. • 

i = l " 

5a. The approximate distance trave4^d is 76 2/3 ^eet. 

b. The approximate wo^k done is ,76 2/3 foot-pounds. 

6a. If n =25. then (-75) (1 + 1/25) + 150 =72. 
if n = 75/ then (-75) (I + I/Z5) + 150 
If n = 300, then (^5)(l + 1/300) + 150*= 7^ 3/^- ' 

b. In terms of the area of the triangle, as n becomes larger, there 

are more rectangles beirig used, with these rectangles forming a 
- ' "closer fit" to the shape of the triangle. 

I f[(J - l)(15/ii))(I5/n) = I ((-2/3)(i -;l)(15/n) + 10)(l5/n) 
1 = 1 . . ^ i = I ^ ' ^ 

= 75 + 75/n. 

b. If n » 25, then 78. 
' If n = 75, then 76^ 

If n = 300, then 75 I /^. . ^ , . 

c. nm(75 + 75/n) = 75, the same numerical value that was obtained 

in (G). , • ^ ' ^ ^ 

8a. Since the midpoint of the ith^ subinterval is given by 

[(i - I)(l5/n) + (i)(l5/n))/2, or, (i)(l5/n) - (15/2n), so 
n 

I f[(i)(I5/n) - (l5/2n)) (15/n) , which = 75 for all values of n. 
1 = 1 ^ 

b. If n =» 25, then 75^- 

If n = 75, then 75. - ^5. ' 

If n f 300, then 75. 21 



c. Iim(75) = 75, the same numerical value that wa^'obtained in (G), 

9. Since the length of each ^ubtnterval is '(8-l)/n = 7/n, then 

,Ax. « 7/n and the left fendpoint is c. = i + (I - I) (7/n), so that 
n 

I f[I + (i - l)(7/n)l-(7/n) = ^*9(1 - 1/n) + kl. 

■ 1 ■ ■ ■ 

10'. ? If^n = 7, ihen ^*9(I - I/7) + '/^z = 8^*. 

If n =» ^9, then 90*' " " ^ - ^ 

i 

etc. 

U. Consider the graph of f(x) =^^2x + A over* [1,^]. 

^ EXAMPLE ^: If n = 7, then the area of the trapezoid in the graph 
^ is approximately 8^ square units. J| 

EXAMPLES: If n =» 7, then from t = 1 to t = Seconds, 8^ feet 
is the approximate distance traveled by a car whose velocity is 
given by v(t) » 2t + ^ feet/sec. ' 
' 'EXAMPLE 6: • If n =7, then from d = I to d = 8* feet. 8^ foot- 
pounds is the afJ^)roximate work done in lifting a bag of salt 
whose, weight at di&tance d T^given t^y 9(d) « 2d + 



12. 



l(2x'+ <<)dx = -ijm • f fjl s['li?-J)(7/n')](7/air; 

- IN I f(l t f^r >'n7/n)](7/.n') ^ 

I3a» Partition [2,5] into n equal lengfhed^'suBinterval s . and select 
' ^ Cj as the rjght endpo^nt. Then Ax. 3/n and c. « 2 + l(3/n), 
for I » 1 , 2^ 3, . . . , n. A Riemann Sum for f(x) « 3x +^5^ 
over [Z,5] is* . ' 

n - 

. ; I f[2 + (i) (3/n)] (3/n) » {'27/2)(l + |/n) + 33. ' 



f (3x + 5)dx - I'im [(27/2) (1 + l/n*) + 33] - 93/2\ 
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EXAMPLE k: 93/2 sq. units is the area of the region bound by 
the lines x-2,x«5, y=0 and y « 3x + 5. 
EXAMPLE 5: 93/2 feet is the distance traveled by a car wh'ose 
velocity is v(t) « 3t + 5 feet/second, when 2 < f< 5. 
EXAMPLE 6: 93/2 foot-poUnds is the work..done in lifting a bag 
of salt^from = 2 feet to d 5 feet, where the weight of the 
salt is 9(d) « 3d + 5. , 

Using the right endpoint of , each equal 1 engthed, subT nterval , 

n ^ ' • 

(5x + I8)dx = Urn I f(-2 + (i)(6/n))(6/n) 
*^-2 n-x» i = l 

• = lim(90)(l + Vn) + = 138. 

\ 

U%ing the right endpoint of each equal 1 engthed- su1>interval , 

2 " 
(x + iix + 5)dx = lim 'J f[-l + (i)(5/n)](5/n), 

= nm[{125/6>(r+l/n)(2+l/n) + 25(l + l/n) + 10] 
n-«» 

• * = 230/3. 
0, as Ax. = 0 for all values of i. 
75/by (G)^ * 

Jn ^(-2/3)t +J0]dt - lim I [C-2/3)(i)(10/n) + 10](lO/n), 
' = linT[(-100/3)(l + 1/n) + 100] 

9' » 200/3. . ' 

EXAMPLHjj : 200/3 is the area of the region bound by the lines 
X - 0, 10, y - 0, and y (-2/3)x + 10. ' 
EXAMPLE 5; 200/3 is the distance traveled by a car whose 
velocity is v(t) » (■;2/3)t + 10 feet/second when 0 < t < 10. 
EXAMPLE 6: 200/3 is the work done In lifting^a bag of s^ft from 
d = 0 to d " 10 and where the weight of the'salt is given by 
9(d) - (-2/3)d + 10. , ^ 

2Z 
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X* 17. /X 

F(x) « ((-2/3)t + 10]dt 



- lim I [(-2/3) (i - '1)(x/n) + 10] (x/n) 

V 

= lim ((-x^/3)(l + i/n) + (2x^/3n) + lOx] 



I8d. 



(-X /3) + lOx, which is the same as (P). 




ength (-2/3)x + 10 



19. 



20. 
21. 



G(x) = (f/2)(x)[I0 + (-2/3)x + 10], or, G(x) = (-x^/3) + lOx, 
which is the same as (P). ^ * 

In (H), let f(t) =* (-2/3)t + 10 be the velocity of a car over 

(0,x]. Since the velocity is the derivative of a position 

function, say fUY = (-x^/3).+ lOx, where F'(x) « (-2/3)x + 10 

which a f (x) . . ^ ' 

I 

a. , b., and c F(8) - F(2) ^ kO. 

(K). rl5 * , - 

((-2/3)x + 10]dx = F(15) - F(0), where .F(x) = (-xV3) + lOx, 



75. 



(l^a.) (k 



I (5x + l8)dx = F(k) - F(-2), where F^x) » (5x^/2) + iBx, 



138. 



(l^Jt).) 

(x^ + i»x + 5)dx » F{k) ' F(-l), where F(x) = x^/3 + 2x^ + 5x, 



- - ^ ^ ^ - 230/3. 
22a. 63; b. 63; c. 63; d. (27/2)' 



e. 677/^. 
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X 



2 

23. The function f(x) = 1/x is not continuous over [-1,6]. 

2k. a. -62 1/2. 

2k 1/8. 

c. Old yoo heed the WARNING in Exercise 23 above? 

d. 8p/9. 

2 



25. 



la. 

b. 
c . 



Area = f (6 - x - x^)dx = J25/6 square units. 

26. Distance = (t^ - 3t + 2)dt + |^ -(t^ - 3t + 2)dt 
+ ^ (t^ - 3t + 2)dt = 11/6. ^ 



27. Since slope = (1^*^* - 72)/(0 - l8) = and the vertical 
intercept is 1^*^*, then q(x) = -kx + 1^^. 

Work = I \kk)dx = 19k!* ft-lbs. 



r 



8. MODEL EXAM 



Set upland simplify a Riemann Sum to approximate the 
area of the region bound by the lines x'= 2, x = 5, 
y « = 0, and y 2x + 3. 

In the simplified Riemartn Sum in la. above, let n-^«>. 

Set up the definite integral that will yield the area 
of the r\3gion in la. abo;t,*e. ' 

d. Evaluate the definite integral in Ic. above by using 
the Fundamental Theorem of Calculus. 

2. Suppose that a particle travels along a straight line 
and its v^elocity is given by v(t) = t^ + 8t + 17 
feet/sec. Find the distance it travels over [1,5] by 
(a) using the limit of a Riebann Sum, and (b) by using 
the Fundamental Theorem of Calculus. 

3. Consider the following limit of a Riemann Sum: 
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lim I ({[1 + (i)(6/n)]^ + 5(1 + (i)(6/n)] - 2}(6/n)) 

Write this as a definite inte'gral and evaluate it, by 
using the. Fundamental Theorem of Calculus. 

Evaluate each of the following i^y using the Funda- 
mental Theorem oL Calculus: 



. (x^ - 3x + l)dx 



-2 




+ x)dx . 




(x + x)dx. 



According to Hooke's Law, the force F required, to 
stretch a spring x" units beyond its natural length*' is 
^(x) = kx, where k^is called the "modulus of the 
spring." Suppose that it takes a 2-pound force to 
stretch a spring from 15 inches (its natural length) 
to 20 inches, so that F(5) = k(5) = 2, or k = 2/5. 
With the same spring, what is the work required to 
stretch the spring from 15 inches to 21 inches?* 

9. 'ANSWERS TO MODEL EXAM 



I {2(2 + (i)(3/n)] + 3}(j/n)» 
1 = 1 

30. 

f5 



I (2x + 3)dx. 



30. 



Ifm I {[I + (i)(Vn)]2 + 8{I + (i)(i/n)] + 17}(Vn) - 616/3 . 
n-»<» is) . # , 



\ 26; 



30 





* b. 




• 
• 


3. 




f 








b. 


8 2V15. 




c. 


f(x) = X 




5. 


1(2/5) 



+ 17)dt = 616/3 feet. 



ERIC • 



27 



31 



umap 



UNIT 426 



MOD' 'IJPS AMD MONfXJRAPHS IN bmiKRORAJPtTATS 
MATHJilAT: :r> and its APPUGATIONS PKOfJECT 



ATMOSPrtERlC PRESSURE IN RELATION TO 



HEIGHT AND TEMPERATURE 



by Arnold J. Ir)$el 



column of. 
atfiiosphere 



X * 



has weight P(x) 



has weight (x) 



ft) 



ERIC 



APPLICATIONS OF CALCULUS TO ATMOSPHERIC PRESSURE 



edC/ umap / SSch^pel st /'newton. mass 02160 



ATMOSPHERIC PRHSSUREJN RELATION 
TO HEIGHT AND TEMPERATURE 



by 

Arnold J. Insel 
Department of Mathematics 
Illinois State University 
Normal, Illinois 6X^61 

\ ■ ■ 

^ TABLE OF CONTExNTS 

1 

INTRODUCTION 

THE SIMPLIFIED MODCL ' . 

2."1 Derivation of the Formula 

2.2 An Example ^ 

THE MORE COMPLICATED MODLL 

3^1 Derivation of the Formula 

5:2' An Example^ T Varies Linearly with x 

3.3 An Application to Meteorology .' . , . 

ANSWERS TO EXERCISES 



Intermodular Descviption Sheet :. VMAP Unit 426 

Title ; ATMOSPHERIC PRESSURE IN •RELATION TO ' 
HEIGHT AND TEMPERATURE " * 

Author ; Arnold J. Insel 

Department of Mathematics 
Illinois State University 
Normal, Illinois 61761 

Review Stage/Date : III 4/16/80 

Classification : APPL CALC/AT^OSHLRIC PRESSURE 

Suggested Support Material . 

Access lt;o a calculator is important . 

Prerequisite Skills : 

1. Knowledge of the calculus as-applied to logarithms and * 
exponentials. , " t % 

2. An acquaintance with high school chemistry is useful. 

Output Skills : 

1. To understand the application of* calculus to construe^ a 
mathematical model of the atmosphere. 




© 1980 EDC/Project UMAP 
All rights reserved. 



r . A • • 

MODITLES.AND MONOGRAPHS IN UNPERGRADUATE 
MATHEMATICS^ AND ITS APPLICATIONS PROJECT (UMAP) ' 

The goal of UMAP is to develop, through a community of users 
and developers, a system of instructional modules m undergraduate 
mathematics 3nd its applications which may be »«fsed to supplement 
existing courses and from which complete courses may eventually \5e 
built. 

The Project is guided by a National Steering Committee of 
mathematicians, scientists, and educators. UMAP is funded by a 
grant from the National Science Foundation to Education Devel#opmont 
Center> Inc., a publicly supported, nonprofit corporation engaged 
in educational research in the U.S. and abroad. 



PROJECT STAFF 

Ross L. Finney 
Solomon Garfunkel 
Felicia DeMay 
Barbara Kelczewski 
Paula M. Santillo 
Donna DiDuci^ 
Janet Webber 
Zachary Zevitas 

NATIONAL STEERING COMMITTE 

W.T. Martin (Chair) 
Steven J. Brams 
Llayron Clarkson 
Ernest J. Henley 
William Hogan 
Donald A. Larson 
William F. Lucas 
R . Duncan Luce 
George M. Miller 
Frederick Mosteller 
Walter E. Sears 
Ge^Jrge Springer 
^miold A. Strassenburg 
Alfre<f B. Willcox 



Director ' 

Consortium Director 

Associate Director 

Coordinator for Materials Production 

Assistant to the Directors* 

Project Secretary 

Word Processor 

Staff Assistant 



M.I.T. 

New York University 

TexarS. Southern University 

University of Houston 

Harvard University 

SUNY at Buffalo ^ 

Cornell University 

Harvard University 

Nassau Comnflinity College 

HarvarcliJniversity 

University of Michigan Press 

Indiana University 

SUNY at Stony. Brook 

Mathematical Association of America 



The Project wodld like to thank Donald F. sfiriner of 
Frostburg State College, Frostburg, Maryland, and Norman V. Schmidt 
of North Iowa Area Coiranunity College', Mason City, Iowa, for their 
reviews^- and all others who assisted in the production ^£ this 
unit. ' ^ 

This material was prepared with the -partial support o'f 
National Science Foundation Grant No. SED76-1961S A02. Recommenda- 
tions expresses are those of the author and do not necessarily 
Reflect the views of the NSF or the copyright holder. 

' 



OCT 



ATMOSPHERIC PRESSURE IN RELATION TO , ^ " 

, ' nmom and temperature 

' -i^'. 1 « iKTRQDUCTION ^ 

'Atjriospheric pressure ort the^'^arth^s surface is due 
r- to the w^>,£ht of the atmosphere above. Imagine a verti- 
r'cal CO iaimn* whose cross section is an inch square and 
n^hich extends lapwards from the Earth's surface witfiout 
bound. (See Figure 1.) ' ' 




r 



er|c ^ ; 3b 



Figure 1^. . * , 

T-he weight of this column in pounds (as weighed in a 
vacuum) is numerically equal to the atmospheric pressure 
in pounds per squ^ire inch at the surface. This is so 
since the pressure at the surface is the force per unit j 
area due to the weight of the atmosphere. More general- 
ly, at^ a height x 'abave the Earth's surface the atmos- 
pheric pressure is numerically equal to the weight of ' 
that portion of the air column above height x. 

We make use of this simple observation ak)ng with 
certain well known facts' about ideal gasses to create '* 
two mathematical models of atmospheric pressure.. The 
first model is somewhat simplified since it does not takje 
temperature variation with altitude into account.^ its 
introduction^ser.\^es the purpose qf preparing t'he way for 
the study of the more complicated second model which does 
take temperature variation into account.. We apply this 
second model to study, condi tions under^-which the.atmo's-. ' * * 
phere is unstable. 



2. THE ^SIMPLIFIED MODEL 



. I Derivation of the Formula 



Consider the column of atmosphere as described in 
the introduction. Let P^ denote the weight of this col- 
umn m pounds. For any,.x > 0 let P(x) denote the atmos- 
pheric pressure in pounds per square inch, x inc'ffes above 
the Earth's surface. Likewise for any x*> 0 let w(x) de- 
note the weight of that portion of the column, in pounds, 
from the Earth's surface to the height of x inches above 
the surface. Figure 2 illustrates the relationship be- 
twe^Jn P(x) and w(x) . . 



has weight P(x)^ 
has weight w(x) 



Figu^2. 

The fol-16wln'g equations are now clear. 
P(0) = P^'' . 



and 
(1) 



>(x) = P^ -Jw,(x) jfor any X >, 0. 



Observe that as^ functions of x, P is decreasing and 
w. is increasing. Let us suppose that the " functions P and 
w are each continuously ^differentiable . 

If we were to weigh samples of air at>various altir 
tudes we would normally find that for a fixed volume, 
samples j)f o^ir., taken at low altitudes are heavier than 
samples ^ttken at higher altitudes. For any x >" 0 let 
P(t)c) denote the weight of a cubic inch of air in the 
cdlumn ^t height x inches above the Earth's surface. We 



lumn at 
all /a] 



shall ^11 p{x) the denaiUj of the atmosphere at height x. 



Ordinarily p is a decreasing function in x. Let 
^^o * tlie density of air at the surface of the 

coiumn. ^ 

Next we assert that for any x >_ 0 ' 

(2) ' w(x) = w' (X) = p(x). 

To see this, first observe that for any^ x > ' 0 and 
Ax > 0, w(x + Ax) - w(x) is the weight in pounds of the 
air in the coliilnn from'height x to height x + Ax. So 
the quotient 

w(x + Ax) - w(x) 
Ax 

is the average density of the air ^n pounds per cubic 

s 

inch in that portion of the column'^ ' Thus 

w»rx1 - w(x + Ax) - w(x) 

represents the density of the atmosphere at' height x, 
thus justifying Equation (2). 

Next, differentiating both sides of Equation (l) and 
applying Equation (2) we obtain *' 

(3)' ' P» (x) = ^-p(x) for any x. 

We now introduce two assumptions used to construct 
the simplified model: 

(a) the chemical composition of the atmosphere is 
uniform and independent of the height. (The 
ratios of the various gasses making up the 
atmosphere are independent of height.) 

(b) . Tyhe temperature of the atmosphere is indepen- 

dent of height*, ' ' 

If w^ apply these assumptions along with the ass^ump- 
tion that tJje atmosphere is^ an ideal gas we may invoke a 
variant of Boyle»s law which states tha.t the density of 

■ • \ I 
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a gas is proportional to its pressure. So for any x > 0 



' 0 
or 
(4) 



POO = P(x) 



P(x) = ^ P(x) 
o 



Combining Equations' (3) and (4) yeilds the equation 
whose solution is evidently 



P*(x)^= P(x) 
o 



(5) P(x) = P^exp 



— p — for any x >^ 0 . 
- o ~i 



Equation (5) iis a formula relating atmospheric pres- 
sure with height. For the sake of practicality we modify 
.(5) so that X is in unit^of feet rather than inches.^ 
Thus (5) yeilds 



(61 ^' P(x) =-P^exp 



2.2 An Example 



--12p^x-l 



for X in units of feet 



Let us now apply Formula (6) based on the assumptions 
of Section 2.1. * ' 

Given the assumptions of this sectipn and that • * 
Po = 14.7 Ibs./sq. inch and p^ = 4.34 x lo'^ Ibs./cu. 
inch, find the atmospheric pressure at 20 ,000 feet $ibove , 
the Earth's surface. By (6) . ^ 

P( 20,000) = (14 . 7) exp [ j Cl2) (4. 54xio'^ ) ^ 

= 7.24 Ibs./sq.. inch. 
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* Exercise 2.1 

/ . , . 

\ Given the assumptions' of this section and the values of P and 
Pq as in the exawple above, at what height is the atmospheric pres- 
sure 1/2 of its value at the surface?- 

Exercise 2.2 ' ^ - * • * 

Given the assuiaptions of this section, suppose that at th^* 
Earth's surface the atmospheric pressure, is 15.00 pounds per square 
inch' while at a height of 1,000 feet above the surface the atmospher- 
ic^ pressure is 14.47 pounds per square inch. Find p^. 

Exercise 2,3 , 

Given* the asSojqptions of this section and the values of and in 
\ Exajr5)le^f'.2, ^^assume. that the Earth»s Radius is 4,000 miles. What is 
the tpiCialyweight of the Earth's atmosphere in^ pounds? 



.5. THE MORE COMPLICATED MODEL 



3?l>*g»^Derivation of the Formula 

Let us now- delete the assumptions thdt atmospheric 
temperature is independent of height. 

We shall measure temperature on the Absolute sc'aie 
of the Kelvin, system. Recall that the Kelvin and the 
Celsius .systems are related since an increase in temper- 
ature of one degree K (Kelvin) is identical to an«ii\- ' 
crease in temperature of one degj^e Celsius, and each 
correspond to a temperature increase of 1.8^ degrees on 
the Fahrenheit scale. However, under the Kelvin system, 
0**^ K is absolute zero, water freezes at 273.1** K, and 
water boils af 373.1** K.^ To convert from the Celsius to 
the Kelvin system, simply add 273.1. 

The relationsfiip between -temperature on* the PaJiren-- 
heit and Kelvin scales are given by the equations ] 



F =^ 1 .8K; 459.58 



and 



K = (F + 459.58)7(1 .8) * 

where F and K are the temperafures on the Fahrenheit and 
, the Kolvin scale respectively. 

FOr any X > 5 let T(x). denote the temperature in 
degrees Kelvin of the atmosphere at height x in inches 
above the'Eartji^s surface. Let^To = T(0) , the tempera- 
ture at the surface. In Section 2 we^applied Boyle»s law 
of ideal gasses to^ obtain Equation (4), In this section 
we complicate the model by superimposing a variation of 
Charleses law of ideal gasses which states that the gas 
density varies inversely with temperature when measured 
on an absolute scale, such' as the Kelvin scale. Conse- 
quently we must introduce the factor Tq/T(x) to the right 
hand side of Equation (4) to obtain 



(7) 



P(x) 



PoToP(x) 
PoT(x) 



Hence 



PITT " PqW 



Hence for any x > 0 we have 



P^(t) 

PltT 



dt 



In 



P(x) 

Pn 



o Q 



T(tT 



dt 



dt 



and^^hen^ce 



?{x) = Poexp 



^0 



JO 



where x is the height in ^nches above the Earth *s surface 
Finally, let us adjust this formula, as in Section 2^ so 
that X is in units of feet. Suppose that the function T 
gives the temperaturje at hei^ght x, where x is in units of 
feet rather than inches. Then the expression T(x/12) 
gives the 'temperatvire at height x, where x is in units of 
inches. Hence at x feet' above the ground 

fl2x 



P(x) = P^exp 



■P T 



0 



1 

T(t/123 



dt 



So we must simplify the expres,sion 
'l2x 



T(t)l2) 



dt 



Introducing the substitution u 
'l2x 



t/1'2 we have 



1 

T(t/12) 



0 



dt = 12 



= 12 



TTtJ 



du 



dt 



Thus we finally arrive at 



(8) 



P(x) = P^exp 



•12p T 

o 



0 



1 

TltT 



dt 



where x is in units of feet and Tfx) is the temperature 
at a height of x feet above the Earth*s surface. 

5.2 An Example, T Varies Linearly with x 

General aviation pilots use the rule: the tempera- 
ture of the atmosphere Mecreaaes linearly with height at 
a rate of 2** C (Celsius) per 1,000 feet of altitude.. 
Since the difference of a degree Celsius is identical to 
the difference <ff a degree Kelvin we may translate thi's 
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rule into the formula 



T(x) = - (2/1000)x 

wh^re X is the height in units of feet above the* Earth's 
surface. Combining this formula with Equation (8) we 
have ^ 



0 



TTtl 



dt = 



0 



1^- (2/1000)t 



dt 



(500) ln(l - x/(500Tq)) , 



and so by (8) we have 

P(x) = PqCxp 
or equivalently 

(9) , ' P(x) = Pod - x/(500To)) 



12p T 

-^^{SOO) In (1 - x/(500T^)) 



6000P T 

00 



Example : ^ 

Assuming Equation (9) with P^ = 14.7 Ibs./sq. inch, 
= 4.34 X 10"^ Ibs./cu. inch and = 293° K 

(a) What is the atmospheric pressure at 20,000 
feet above the Earth's surface? 

(b) -Jtfc:W^at;--^£<#ght is the a tmosp her i c f^ re s sure " 

half o f rhr_ pr rTTffrT n t the. ■ iirtn r rf 

Solution of (a) : 



P = 14.7 



1 



20-,'000 
(500)(2d3) 



(600O)(4.54 X 10'^)(293) 
14.7 



= 6.86 Ibs./sq. inch.^ * 
Comparing this answer to the answer o^ 7.24 Ibs./sq. inch 



i3 



8 



of the problem as treated in Section 2.2 we see tha-t tjiere 
is a difference of only about S.S%. This i's not surpris- 
ing if we observe that the first order Taylor approxima- 
tion of Equation (9) is indepej^dent of T^ . 

Solution of (b) : , — 

^ Assuming P(x)/Pq = 1/2 Equation (9) becomes 

')(293) 



= 1 - 



[500X293) 



(6000)(4.34 X 10" 
TO 



So 



i 

T7T9 



1 - 



= 1 



(500) (293) 



(S00)(29S) 



5.19 (computing to three 
significant digits) 





/ 




1 (5.15) 


/ 


X = (500) (293) 


1 - 


2 



= 18,300 feet (to three significant digits) 
3.3 An Application to Meteorology 

! 

Let p(x) denote the density of the atmosphere at' 
height X in feet. We would normal ly expect that p de- 
creases as X increases. Under these conditions we shall 
say that the atmosphere is stable,^ Otherwise we shall 
call the' atmosphere unstable. In «n unstable atmpsph^re 
"a given volume of air above wo^ld weigh at least las much 
as an equal volume of air below. Under these circum- 4,* 
stances there would be a vertical motion of air causing 
winds and down-draughts. Mathematically the air is un- 
stable if 



(10) 



P*(x) > 0. 



Let us study the conditions of instability under the 
assumptions and results of, Section 3.1. Combining ^ 
Equations (7) and (8) we have 



p(x) . ^ exp 



o o 



lo ntj 



dt: 



where x is in units of feet and T(t) is the temperature 
(measured on the Kelvin scale) t feet above the surface. 
Thus by Equation (10) the atmosphere is unstable if 



P.;(x) 



(T(x))2 



T' (X) 



x) + 



o o 



exp 



This inequality reduces to 

-12p T 



T' (x) < 



Q Q 



■12 r 

O 0 



TJJJ 



dt 



> 0. 



Exampl e 

Assume that T drops linearly with height, Tq =' 
293*! Ks pQ = 4.34 x iq-^ ib$./cu. inch, and ?^ ^ 14. J 
Ibs./sq. inch. Fi,nd the maximum^ temperature at 1000 feet 
above the suffate of the Earth so that the atmosphere is 
unstable. The assumption of linearity requires that. 



T(x) 



kx 



■for some positive constant k. 
have 

--12p T 
-k < 



Thus by inequality (11) we 



or 



0.0104 £ k (approximating , to three 
.significant digits) 



To achieve the maximum temperature we require that k be 
as small as pcrssible, i.e. k = 0.0104. Thus 

T(x) = 293 ' (0.0104)x. 

So for X ^ 1,000 

T(x) = T(1,000) = 293 - 10.4 = 282.6. 

Notice that there is a temperature drop of 10.4 degrees 
Celsius per b,000 feet. This corresponds .to a dVop of 
approximately 18.7 degrees Fahrenheit per 1^000 feet. 

One tinal note. The second model is still over-sim- 
plified since it does not consider the' possibility of 
variation of atmospheric composition with altitude. In 
particular, tfie atmosphere* may vary in altitude with rec- 
spect to the amount of water vapor it contains. The 
atmosphere may al«6o contain such po^llutants as smoke and 
smog. These all contribute to its density, and hence to 
Its pressure. 

Exeiyise 3.1 ^ ' 

Assuming Equation (9) , = 14 . 7, = 4. 34 x 10-^ ^ an^l Tq = - 
300* K, what is the atmospheric pressure at 10,000 feet and at 
20,000 feet? ^ , ' • 

Exercise' 3. 2 * ' i 

Assume the model of this section and that temperature decreases 

linearly with height. ' Suppose Pq = 14.7, and T^ = 283. Suppose, 

in addition, that at 10,000 fqet T = 293 and P =.10.2, Find Pq. 

>* ♦ ♦ 
Exercise 3.3 

Assuming T is a constant (and therefore T(x) = T^. for all x) 
show that Equation ,(8) reduces to Equation (6) . 

Exercise 3.4 



A|si^ng that T decreases linearly with height and assuming the y 
raod^l of^this section, generalize tfie formula of Equation (9), giveri 
that decreases by k degrees Celsius per 1,000 feet. 



4. ANSWERS fo'^ELECTE.D EXERCIS£S 

19,600 feet. ' , 

= 4.497 X 10"^ 
1.19 X ' ' 

10.2 lbs. per square inch at 10,000 feet. 
6.87 lbs. per" square inch at 20,000 feet. 

p = 4.40 X 10"^ 

0 



P(x) = P, 



1 - 



kx 
lOOOT^ 



12000P T 

o_o 

P k 

0 
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J ^Not enough detail to understand "the unit 
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^Appropriate amount of detail • .'^ - 

JJnit was occasionally too detailed, but this was not distracting 
Too much detail; I was often distracted 



2. How helpful were the problem answirs ? . , ^ 



Sample solutions were too brief; I could not do the inte^JSiate steps 
Sufficient ^Information was given to solve the problems - 
Sample solutions were too detai3^ed; I didn't need ther^ 



3. Except for fulfilling the prerequisites, how much did you use other sources (for 
example, instructor, friends, or other books) in order t6 understand the unifc? 



A Lot ^Somewhat A Little Not at all 
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« as many as apply,) , 
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Paragraph headings ^ , - 
^Examples * 
^Special Assistance Supplement (if. present) 
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Were any of the following paVts of the unit particularly helpful? (Check Is many 
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^Statement of skills and concepts (objectives) 



_^Examples 
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Please describe anything in the unit that you did not particularly like. 
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^ 1. INTRODUCTION 

In this module we consider some properties and illus- 
trations of space curves that lie on-'given surfaces. We 
explore the following' questions: 

1) How IS the instantaneous rate of climb along a 
curve on the surf/ce related to the equation of 
the surface? ^ 

2) Is It possib^le to find a path on the surface 
for whi;rh the rat^ of climb is optimized at 
each ppint along the way? 

— \ 

As ph)^sical exai4)le?, consider the^fpHowing situa- 
tions : ^ ' 

' E^cample 0-1 . A party of mountaineers can ascend to 
a summit by a leisurely proce&s of "swi tchbacking" . A 
few hardier members of the group wish to test their met- 
tle by^starting from the same location as the others, but 
reaching the summit by .a path that is always the steepest 
possible. Certainly s.uch a course, if it qan be found, 
will be much more strenuous! If the topography of the 
motintain is known, the principles of the gradient may be 
used to ch^rt such a course. \ 

Example An oil tanker Jias- met with disaster at 

night,, on a calm sea, and is left without radio communi- 
catlo^. A rescue vessel that is able to monitor continu- 
ously the con^:entration .C of the spreading oil slick 
tri.es to locate the tanker by moVing in the direction of 
greatest in&Yease opthe concentration* What is its path? 

In the situation of^Example 0-2, we are thinking of 
the concje'ntrAtion of oii^:, ^C, as the dependent variable, 
and we^^have C = f (x ,y) Vhich<;;is of the form z = f(x,y),' . 
(a standard djesignation for a 3-dimensi'onal surface). ' : 



2. 



PRELIMINARIES 



'Let us begin by considei;ing\what is meant by ''a 
curve lying on a si^rface". * Suppose S is a surface in 
three-dimensional space E3 defined by 

'/^^ f(X,y) for all^(x,y) in D , 
where D is' a set, in the xy-plane. We ^assume that S is a 
"smooth" surface; i.e., the partial derivatives 3f/3x and 
9f/9y are continuous at each point in D/ 

Let be a durve that is defined by the parametric 
equations ^ - * 



(A) 



^x = x(tr 

y = y(t)' 

2 = 2(t) 



where t varies over an interval I, on the real line. If 
x(t) and y(t5 lie in D for ca'ch t in I, and if z(t) sat- 
isfies the equation z(t) = f [x ( t) , y ( t) ] , then we say that 
y lies dn the surface S. See Figure 1. 




Figure 1. 



We also stipulate that Y be a dif ferentiable curve; 
i.e., y will possess ai:c length for all finite intervals 
(tj^t^), in I. This is equiA^ent to the requirement 
that the funct ion rx( t ) 1 ^ + ly(t)]^ + (z(t)]^ be contin- 
uous for all t inpt,,t,) and the arclength s of Y from 
tj to tj^ be given by 

s = /(x(t)]' + Iy(t)]' + I=(t)]^ dt . 
< t 



The symbols x(t), y(t), 2{t) Indicate, respectively, 
dx/dt, dy/dt, dz/dt. 

In applications, D may be a bounded or unbounded do^' 
main in the xy-plane and I may be a finite or infinite 
int^erval. • 

Example 1 . Consider the surface S givert by 
(1) z = + xy + Zy"* , 



or 



Note that 



f(x,y) = x^ + xy + Zy"* 



3f - -„ . 0,.3 



3y = ^ * 8y 



One specific curve y lying on S is given by 



Y: 



X = x(y) = t 
y = y(t) = t^ t in 11,2] 

z = z(t) = t^ + t^ + 2t® 



where we mav verify that z(t) « [ x{t)]2 + x(t)y(t) + 
2Iy{t)l** by i sub^stituting into Equation (1). The arc- 
length s'oythis curve is given by • \ 



(The value of s can be found by numerical integration; 
we used a Romberg proce4ure to obtai^n s - 520.0302 .) 



Exercise 1 : Let S be the surface, defined by 
2 = f(x,y) = x^ + 2yS for all (x,y)'. 

a) Sketch the portion of S that lies over the Cirst octant. 

b) ' Find 3f/ax, 3f/8y .\ 



c) Verify that the curve Y*. 



x(t) = t cos t 
y(t) = ^= sin t 



2(t). 



d) Find the arclength of the curve y *for t in' [0,2tt], 

(Hint: You should see that Y is a helix winding up around the 
elliptical paraboloid 2 = x^ + 2y^.) 



5. ^DEFINITION QF THE GRADIENT 



The partial derivatives of f(x,^) determine aij im- 
portant vector field in the xy-plane.; 

•The Gradient Vf(xQ,yQ) of the surface S defined by 
the function z = f(Xjy) at any point (x^,y^) in D is the 
vector 



Vf(Xp,yo) ^ H i - |~ 3 



9y 



where i and j are unit vectors along ^the positive x-axis 
and the positive y-axis, respectively, and the notation 



\ 



(Xo>yo) 

indicates that we substitute x^ for x and.yo for y in the 
expressions 3f/9x and 3f/3y. 
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Example 2 . To find the Gradient of the function fn 
Example 1 ,at the point (x^,y^) = (1,2), we find the par- 
tial derivlatives . \ « 

II = 2x . y , II = X . 8y^ , ■ 
form the general expression 

^f(Xo,yo) = (2x + y)i + (x + 8vMj 

(1,2] 

then substitute 1 for x and 2 for y tof find 

. Vf(l,2) = 4i + 65j 

Note: The Gradient S depends only on f(x,y) and on 
(Xo,yo) but not on z^, where = f(Xo,yJ. This fact, 
though readily apparent, deserves emphasis! Even though 
S is d-dimensionaly the gradient of S is a 2-dimensional 
. vector. Often the dradient VfCx,y) is confused, with the 
Normal to S at the point i^^^^V^yZ^)^ The Normal 
^^S^^o ''^o) ^» however, is given by 

which is a vector with three components, while the Gra- 
dient has only two. 

4. -'^-FINDING THE TANGENT VECTOR TO A CURVE ON ^ 
^ ^ A SURFACE USING-THE GRADIENT 

Suppose that a surface S is defined by z = f(x,y), ' 
and that y is a curve on S. Then recalling the paramet- 
ric form (A) for 7 we may defjne T(tQ), the tangent vec- 
tor to Y at t = tg by the equation 

f(tj - i(t)i + y(t)j + z\x)\i i ' 
" ^ ^ ft = 

z(t) |nay be found directly from the parametric forms, 
but since we have a chain 



f(x,y) : 



X = x(t) 

y = y(t) 



er|c • 



we may also use the .appropr j.ate chain rule, as follows 



(BO 2(to) = II [x(t),y(t)]i(ty + |i (x(t),y(t)]y(t3| 



t = t„ 



or,*what is the same thing 
(x,yj 

x=x(tj 

y=y(tj 



z(t. 



ax 



x=x(tj 

y=y(tj 



y(to) 



The chain rule is also commonly written as 



dz _ 3f dx ^ 3f d 



HT ~ ax 3T 9y 

see Thomas 5 Finney, p; loi or Greenspan 5 Benney, p, 496 
An example should help t9 clarify the two distinct ways 
of finding z (t) . 

Example 3 . From Example 1, we know that the curve 
X = x(t) = t 

y = y(t) = t" t in [1,2] 

z = z(t) = t^ + t^ + 2t^ ' 

liesonthesurface 

z = f(x,y) = x^ + xy + 2y** . 
From the definition of y we have 

z(t„) = 2t„ ^ 3t„^ . 16t/ ; 
SO that tangent vector to y , at t^ is given by 

I T(tJ = Ti . 2tJ . (2t^ ^ 3t^2 . 16t/)k 

can also calculate ^{t^) by using the chain rule, a.^f 
follows : 



35f*(x(tJ,y(to 



)} = 2x + y 



(x = to,y = to^) 



= 2t^ + t/ 
0 0 



H (x(tj.y(tj) = X,* 8y' 



from which, using the chain rule (B) , 

Hence, the tangent vector to > at t^ is given by 

T(to) = li + 2tJ + (2to + 3t/ + 16to')k , 
which checks with the above calculation. 

Exercise 2, For the curve 



7: 



' x{t) = t + 1 
y(t) = t - 1 
z{t) = j 



1 



which lies on the surface z = xy, find the tangent vector T(l) in 
two distinct ways. 



The connection between the Gradient and this 
alternate calculation which uses the Chain Rule will 
be explored after we recall* some concepts and proper- 
ties of vectors. 

Let a be the vector a = a^i + a^j + ^3^^- ' ' 

1, The length of the vector a is denoted by ia| and' 
*is given by * 



2. The projection of a on the xy plane wi^ll be 
denoted' by * 

^p ^ * ' Figure 2a. 

3. The dot product of a with vector = b^i + b^j + 
b^k is denoted by a • b° and is given by a ' b = 
a^bj + a^b^ + ajbj. We' also have a S = ia| |Sf 
cos e, where 6 is the angle betweeij a and B; see 
Figure 2b . , 

4. Two non-zero vectors a and b are said to be' 
orthogonal JioT perpendicular.) if a ' S = 0, i.e., 
if the angle G between them is t\/2. From 3, a 
condition for orthogonality 'also becomes*a,b- + 
a^b^ + a^bg = 0, which is satisfied in two dimen- 
sions (note a^ = 0, b^ = 0) 'if and only if we have 
b = cC-a^i + aj), where c is a proportionality con 
stant. * ^ 

So far we haye seen that if the curve y lies on the 
surface (given by z = f(x,y)) and if x = x(t) and 
y = y(t](^re given, then the z-coordinate for y must 
satisfy the""surface requirement" z(t) = f [x(t) ,y (t) ] . 

More impoxtantly we notice that, as in (B) , 

z(t) ■= ^ [x(t)fy(t)]i(t)|+M (x(t),y(t)3y(t.) 
I 

which may be written as the dot product 

I z(t) = Vf(x(t).y(t)) -.TpCt), 

where tp(t) = x(t)i + y(t)j is the projection of the 
tangent vector«f(t) in the xy-plane. 

* Gl 



I 




Thus the rate^of change of the coordinate z(t) de- 
pends, on two vectors, Vf (x(t) ,y(t)) and Tp(t) , both of 
which^saj^e located in the xy-plane. By the definition of 
the dot product we can also write 



^^2(to) = tvf(x(t),y(t))| |Tp(t)| cos<|;(t) 



wher^ii;(tj is the angle between the Gradient and the 
projected vector Tp(t) at t = t^, . We see, then, that 
the sign of zit^) is deterlined by the sign of cosij^Ct) 
and thus that the coordinate 2(t) increases when 
cosii;(t) > 0 and decreases when cosij^^t) < 0. 

And so i(t^) may be found from the two vectors 
in the xy-plane (Figure 3) merely by forming their dot 
product. * 

Of particular significance is the situation where 
'l'(t) = tt/2 for all val^ues^of t.in some interval. In jthat 
case, cosi|;(t) s 0, and z(t) (from above) for all t 

in the interval. Since the rate of change of z(t) is 0, 
z(t) is constant for all t in the interval. This leads 
us to consider only the family of level curves , 

y-" f(x,y) = c, or f[x(t),y(t)] = c 

which are projections of curves y lying on S such that • 
each curve y is parallel to' the xy-p,^yyie.-(all along the 
curve Y, z has' the same value c) . The notation y, will 
be used to refer to a particular level curve of the form 
f(x,y) = c; see Figure 4- 

Thus cos<i^(t) = Q implies: 

At^each point (x ^) on a level curve f(x,y) « Cq / 
^n^o^Yo^ • 'Tp(t J - 0 

where Xq = x(to), yo = /(t^), and T^(t^) is constructed 
from the parametric*equations of the level curve (in the 
case they are available). If parametric equatio^us for 
the level curve f(x,y,) = are not readily available, a 
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, ' ' ' i 




Figure 3 




Figure 4. ^ 

vector proportional to T^Ct^) may be constructed simply 
by noticing that ^fCx^.y^) and ^pCt^) are ^er^pendicular- 
so' that by property 4 above, 

Some examples will clarify these .'concepts . - ' 

Example 4 > One level curve Y^^ for the surface S: 
z = x^ + 4y^ is S = x^ + 4y^ , and it contains the point 
'(x^,yj = (1,1)' See Figure 5, For this level curve, 
we .write some parametric equations with relati'oe .ease: 



G5 




V / 



fx(t) = /5 



cos t 



y(t) = sin t , 



t in (0,2Trl 



Figure Sa. 




Figure 5b. 



13 



and at {x^,y^) =^(i,i), cos t^, =Na//5, sin = 2//5. 
Now 

Vf(x^,yJ = 2x i ^•Sy 1^ = Zi + 8j;* 



sm t^i + — I cos t^j 



■ -.-^S • I i - I (-41 . 5);.. 

And we check that on this level curve, 

^ff^o'^o) '^^'^p(S) = > 8j) • \ (.8i ^ 2j) 



(-16 + 16) = t) 



Example 5 . Next, treat the same surface a5 in Exam 
pie 3, but without parairfetrizing Yj^. The level curve 
5 = + 4y2 contains the poirit (1,1). As before, 

^f(x^,yj •= 2i > S] , 

and by in'spection we construct 

Tp(tJ = c(-8i ^ 2j] , 

which hits the same direction as the ^) in the pre- 

vious exan^ple. ^ 



I 



5 . SUMMARY 



^ ^Let (xQ,yp) be a point on Ithe level curte Y^^: 
f(x,/) = Cp. Then the vector * 

is normal to y at (x ,y ), and the tangent vector to Y. 



at (x,,yj = (x(t;) .yCtj)) is 
or 

A final example will illustrate the use, of the Gra- 
dient m the case where the level curve is not ,easy to 
parame tri ze . 

Kxample 6 . Let z = f(x,y) = x^ + xy + 2v\ A Jevel 
curve Vj. for this surface S is + xy + 2y'* = 4, and 
C^o'^o^ = (A»l) IS a point on y^. As can be seen, is 
not easy to parametrize, but from Example 1 and the above 
si^mmary, we may^write 7f(x^,yg) = *5i + 9j and, by con- 
struction, Tpit^) = (-51 + j). iTrrfTtrtrtalTy , this latter 
could also be obtained by noticing thaT dt^^t differen- 
tiation of the level curve yields 

(2) (2x + y)dx + (x + 8y^)dy = 0 ' 

or, at the point (1,1), 

' 5dx + ^^^^ ^ * 

thus forcing the tangent vector, ' x, 

dx ^ dy c ■ f 

at- ^ * dt J ■ - . ' ■ . 

to have a directi^ft parallel to '3i + J- ' ^ ' 

Exercise 3 . A point P is nxjv^g along the level curve of Example t\ 
At a Certain nwinent its acceleration l^s known to be a = 4i - 2o , and 
it is located at the point (1,1). Resolve the acceleration vector 
into two components, a tangential component and. a component normal 
to the 'path (Hint: If a vector 'l^ is. resolved into tairgcntial and 
normal components T and N, such as a = CjT C2N, then taking the 
dot product vof both sides with T fields *f • a = Cjf • f since 
f • ii i= 0; thus • . , 

" ' .IS 



and in a 'similar manner we fin^ 
/ fJ • a , 



Question : How does Equation (2), p. 20 differ from the 
chain rule (B) on page b for i(t) mentioned earlier in 
this sec t\on? 

Review Question : What is the role of the Gradient in 
constructing the tangent vector f(t) to a curve lying on 
a surface S? - ^ ' 

'Answer: ^i) for an \arbitrary curve v , the planar compo- 
nents x(t) and y(t) would have to be given, from wrfich 
the components x(t) and y(t) are obtained; then the Gra- 
dient may be o^ed (via the chain rule (B)) to determine 
the colniponenH: 2(H^.' 

^y^Tvi for a level curve y^, we may find a vector 
* proportional to T(t) by merely choosing a vector perpen- 
dicular to the Gradient vector. 



6. APPLICATION OF THE GRADIENT TO FINDING 
CURVES OF STEEPEST ASCENT (DESCENT) 

Cocisider again the preliminary example mentioned m 
the introduction. Suppose that f(x,y) m Figure 6 repre- 
sents the concentration of oil, C = f(x,y). Now if the 
rescue shijj^^ is somehow able to monitor the concentration 
of oil at each point (x^jy^) and wishes to go in the di- 
recTion of greatest increase of C (from left to right) 
proceeding from level curve to level curve, we would ex- 
pect 'from Lihe results of Section 4 that it- would contin- 
uously follow the Gj*'adient direction from the point, say 
(Xo,yo), to the right. Actually, this is true, and we 
will provide. a proof^of this before continuing with some 
examples. - ^ 

CD 



Figure 6. 



Keep in mind that we are seeking a path in the 
xy-plane which will lead the rescue ship to the disa.bled 
tanker by the method described in the preceding 'paragraph 
Let us translate our goal into more forona^ mathematical y 
language. The desired* path in the xy-plane could be in[ 
terpreted as the projection Yp(t) onto the xy-plane of k 
special curve Y(t) which lies on the surface C = f{x,y) 
<5r z = f(x,y), (where we see that the roles of C and z 
are interchangeable!). The special ^rof^erty of Y(t) is 
that its projection Yp(t) in the xy-plane, continually 
moves in the direction of greatest increase of C = f(x,y) 
from level curve to level curve* 
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Our method of solution will be. to discover the curve 
Y(t) and then to return attention to its projection y (t) 
in the xy-plane, which is the path we seek. 

The cjurve Y(t) can be described parametrically as 
follows: 

fx = X(t) 

T(t) : ■ y = Y(t) a < t < b 

^= Z(t) 

here, since yit) lies on the surface C = f(x,y) (or z = 
f(x,y)), we have Z(t) = f [X( t) ,Y (t) ] , and the projected ^ 
cjurve Tp(t), a» in Section IV, is*^ 

[x = X(t) 

Tp(t) :l a < t < b 

(y = Y(t) 

As before, the symbol Tp(t) denotes the tangent vector to 
the curve Tp(t), and as in (B) , we also have 

(C) 2(t) = ^f [X(t),Y(t)] • Tp(t) . 

As you are aware, any curve may have several (and 
even infinitely manyl) different parame tri zations ; we 
wish to avoid IrTy complications which may arise *along 
this line by scLnehow ' "normali zing" the curve Yp(t). It ' 
is also -for convenience that we wish to make Tp(t) inde- 
pendent of the parameter t, and thus we introduce the 
arclength s which is independent of the particular pai:a- 
metrization used for y^it) . 

(D) • s = |^Tp(t)idt = |V[X(t)]^ ^ [t(t)]^ dt ' 

0 Q 

or (by the* Fundamental Theorem of Calculus) 

ai = I'^pf^)! = /[x(t)]^ - [Y(t)]^ . 

Note^|hat, regardless of the form of parame tri zation of 
Y(t), we can solve (D) (symbolically!) for t as a func- 
tion of s, say t = t(s). Now, dividing both sides of (C) 

18 

71 , 



by |Tp(t) I , wc obtain 



iii^ = (X(t),Y(t)) • 4^^^^ 



l'Fp(t)[ » |Tp(g| 



or 



dZ 

ar dz 



^ = (X(s),Y(s)] -.11^^15) 



ds 

at 



where 



Uj (s) s 



is (conveniently I ) a unit vector in the direction of 
^p(t), and we have replaced t by the expression t = t(s) 
which was found (symbolically) from (D). 

Once again we recall the rUle^s for a dot product and 
note that 

dZ 



l^f (X(s),Y(s)] I |ut^(s)| cos^(s) 



[^f (X(s) ,Y(s)] j cosiI;(s) . 



What have we achieved here? This expression for dZ/ds 
relates the rate of change of Z with respect to distance 
(arc length) along a curve lying on a surface S and the 
Gradient of that surface, [ 

From this expression, we can sjee that 

(i) dZ/ds IS maximized or minimized by allowing 
^(s) to be 0 or tt , resf^ectively , and 

(ii) max dZ/ds « 1 7f [X ( s) , Y ( s) ] t , 

min dZ/ds = - | ^f (X (s )„ Y^s) ] | ; 

either (i) or (ii) wilt be attained if the direction of 

the unit vector u-r (s) is precisely the same (or oppo- 
P 

sitel) direction as the Gradient; more precisely/ if 
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(s) = . ;f[X(s),Y(s)] 
P j |Vf[X(s),Y(s)].| 

I t • 

Thus^ we see that our assertion in the] first para- 
graph of Section. is indeed valid: the path to follow 
to achieve the greatest local increase (decrease) m 
C = f(x,y) IS precisely a path which takes the direction 
of the Gradient vector (or opposite that direction) at 
any point (X,Y) = (X(s),Y(sJ]: i.e., at each point the 
decision on which direction to move next is made on the 
basis of examining the Gradient. We will refer to this 
path as the "path of s^teepest ascent (descent)**. 

Another way of stating this important result is to ' 
note that the tangent vector uj (s) to the path of steep- 
est ascent (descent) Yp(t) is parallel to the Gradient, 
Vf ()i(s) ,Y(s) ] ; thus the components dX/ds and dY/ds,of the 
tangent vector can be expressed in direct proportion to 
those of the Gradient. 

^ = X 1^ (X(s),YCs)] , 

g = X II [X(s),Y(s)J , 

where *X = X*(s) is a function of proportionality. 

Example 7 . If B /^and A^^^. then t;he level 
curves of the surface S: 2 = f(x) =Vax + By which ' 

is a plane, are straight lines 

Ax + By = C , 

each of whose slopes is m = -A/B. We will construct 
Yp(t) , the curve of steepest ascent in the xV-plane 
through the point (0,0), and Y(t) the curve Of steepest! 
asc^nt^lying in the plane S. We have 

(3) dX . ,. dY . 
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an4 since the curve is to*pass through the point (0,0) we 
may choose to measure arc length s from that j)oint, so 
that an additional condition on the curve of steepest as- 
cent becomes 

i^) x(o) = 0 , y)o)'= 0 

Then from (3) we have 

dY . * 
• JlX " .5 * 

Now integrate and apply condi*tions (4). 
Y(s) = I X(s) 

Hence we see that a set of parametric equations for the 
curve of steepest ascent in the xy-plane is. 



' X(s) ^ s 



B s > 0 

Y(s) = ^ s 



and a set of parametric equations for the curve of steep- 
est ascent Jlying on the surface S is I 



Y(s): 



X(s) 
Y(s) 

2(5) = ("A + ^)s + D 



Y(s) = I s s > 0 V. 



It is now evident tha't Yp(t) (where we replace -the para- 
meter s with the parameter t to facilitate correlation 
with concepts developed earlier) is perpendicular to every 
lever curve in the xy-plane. This is no accident,' an^ 
we leave it to you ^to prove: 

A steepest ascent curve Yp(t) is perpendicular' 
(orthogonal) to any level curve at a point of 
intersj^tion . 

Example S \ Find both steepest ascent cvrves Y(s) 
and Yp(s)^for the surface z = ^Cx^ + y2)>through the 
point (Xo,y, ). ' j 
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- ^x 



dY 



dY Y 



fx = Y 



Upon separating variables, 
dY dX 

and integrating, mY = mx + inC = inCX or Y = CX. Also 
since Y = |y„ when X = Xo ve have Y = (yo/Xo)X.. 
A parametrization for Yp(s) would bef 



Yp(s) 



and 



Y(s) 



[ X(s) ^ s 
' X(s) =. s 



0 ^ s < X 



0 1 s <• X. 



You will notice that "Gradient" curves Yp(s) are straight 
lines whenever the level curves f(x,y) = C are circles 
(see Figure 7); keep this in mind as you proceed to the 
next example. 

Example 9. Find both "Gradient" curves y(s) and 
Yp(s) through the point (X,,yJ when z = f(x,y) -f^ 
h (x^ + uyM , u > 1- 



dX 



dY . 



3? 



= UY 



dY _ uY 

3x " X 



froip which we form the differential equationj 

^ = U ^ ; X(b) = X, , Y(0) = y^ 
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Integrating, we obtain 

InY = plnX + InC 
= InX^ + InC 
= InCX^ 



Y ^.CX'' 

Again inserting the point (x ,v ) we have Y = v (X/x )V. 
Thus paralmetrizations for yis) and Yp(s), respectively, 
would be 



f X(s) 
Y(s) 



0 X, 



0 •< s < x„ 



{ 



Yts): 



X(s) 



y(s) = y„l,^ 



(f^< S < X, 



z(s) = i(s^ . -'yo'i^V] 



In the present case, we see that the steepest ascent 
curves in the xy-plane resemble parabolas passing through 
the origin whenever the level curves of f(x,y) = C are 
ellipses centered at the origin. 



Exercise 4 . Sketch the level curves and the gradient curves Yp(s) 
for the previous example in the cases where p = 2 and p = 3. 

f 

i ■ 

A remark concerning the solution of the differential 
equations is in order. In finding the gradient cui;ve 
Yn(s), you will always be using expressions of the form 



g = f,(x.y). 



3? = fy(^'Y)'. 



aild consequently 
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-ax - f^(x,Y) • .• 

The general problem of solving this differential equation 
for the actual curve ,Y = Y(X) which we have been refer- 
ring to as Yp(s) can 'be quite formidable. In our e/am- 
pies, you will notice, the expression (Ej gave rise to a 
rather simple separable differential equation. In every 
case, this was 'the result of .special choices for f(\,y), 
and a little imaginat ion ' shouj^jd (^vince you that a more 
complicated sur^ce z = f(x,y) will result in a more chal- 
lenging differential equation. Along these lines, we sug- 
gest" that those of^you who are better versed in solvinTg 
first-order differential equations explore the fcxjlowmg: 

^ Problem : Find the Gradient curves for the surface 



where AC - > 0, B > 0, such that Yp(s) passes through 
^ "".^j* y " Vq' ^ou will find that (E) becomes a homo- 
geneous differential equation. A further hint: Perhaps 
an initial rotation of axes to eliminate- the xy term 
would prove convenient! 

We conclude with an additional exercise which re- 
lates our development again to a physical setting. 

Exercise 5 . A tanker located at coordinates (x, ,yj) has capsized^ 
leaving an oil slick floating on the calm surface of the ocean with 
concentration given by the law 

C(x,y) =^C,e-^I^^-^i^'^^^^^y->^'^ . > 0, .>!. 

/ 

A small rescue vessel proceeding from location (x^.y^) moves slov|y 
in the direction of tM^fr^asiM^^ concentration according to the above 
law for a calm sea. What is tl^e equation of the rescue vessels 
path? , ' ♦ 
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7. MODEL EXAM 

A ysugarloafl' mountain has the equation 
2 = H - cx(x' + yy2), a > 0, y > 1. Find the steepest ' 
ascent' curve a mountain climber must take if he or she 
wishes to ascend the mountain by Gradient methodsl The 
climber starts at the point (x^,y^,0). (Follow the indi- 
cated steps to the solution.) 

a) Find the level curves for the surface 

2 = f(x,y),. What type of curves are they? 

b) Find 'the Gradient vector for the surface at the 
point (1,2), and also at an arbitrary point 
(x,y) 

c) Firtd the 'level curve passing through the point 
(1,2). 

d) Find a tangent vector"to the level curve in 
c) at the point (1,2). 

e) Using the result of b), find the differential 
equation for the steepest ascent curve. 

f) Solve the differential equation insertin'g the 
given condition a t X = x^ , Y '= y ^ . 

g) Find parametric expressions for the steepest 
ascent curve y(s) which actually lies on the 
mountain surface . 

h) Choose your own values for a and y , and make " 
sketches of the mountain surface, the level 
curves, and the steepest ascent curve y (s) . 
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8. ANSWERS TO EXERCISES AND MODEL EXAM 



Exercises 



a) The surface is an elliptical parabqloi'd, with 
lowest point at the origin; its sketch i? very 
simi'lar to that of Figure 5a. 

h) . = 2x , . ^ = 4y ^ 



c) 



lx(t))' +^2(v(tj)^ = t^cos^t + Z^irj sin't 
^ = (t^ cos^t + sin't) 

/ - t' = z{t) 

d) Using the formula for arc length given in. 
Section II, and grouping terms,* / 

t=2TT • • ^ 



s = 
t-0 



( * ij^in^t + + ijcosH - t sint co^t +4t' dt 



From- the definition of y and T(to), with. to = 1, we 
have, ' . 



T(l) =' i(t)i ^ y(tlj + htjk 



t = 



or T(l) = li f Ij + 2k. Alternately, we can use the 

chain rule; here f(x,y)'= xy. Since x(l) = 2 and 

XI) - 0, • » ' ^ 

||(x(i),v(i)) f - 0 . 

ly = 0 

and • 

* ilcxdJ.Vd)) = xj^ ^ = 2 . 

j X — *. 

y = 
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from which 
- * 2(1) = 0 x(l) + 2y(l) 

0(1) + 2(1) = 2 

so that 

T(l) = li + Ij + 2k . 

4. From analytic geometry, you should find that the 

level c-urves for both the cases y = 2 and y = 3 are 
m 

ellipses, with the imajor and minor axes maintaining 
a constant proportion to each" other in each'case^ 

The gradient curves Yp(s) are parabolas for the case 
* y = 2 and cubic cur/ves for the cas^ y = 3. Sketch 
the corresponding families* of curves on the, same 
graph, observing the ortho'gonali ty property! ( 

*dX C^(-2a(x-Xj)) 
ds " M 

C J-.2ay(y-^yj>) , 
ds " ^ 

we have used M to^ represent 

^-a[Cx-x^)2 + u(y-yj)^] . 

dY ^ ^^y^y'i^ , 

3X " x-x J ' 



5. 



where 



Then 



81 



hj For simplicity, let us choose a = 1, u = *4 . 

.2 = f(x,y) =:H - (x^ + 4yM, an elliptical para- 
boloid. The level curves are H - C = ^. 4y^ , 
0^ 1 C- 1 )f; a family of ellipses. 

The steepest ascent curve Yn(s) is ' 



^See Figures ?^a, '8b, and 8c. 
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separating variables we find 
M dX 



dY 



x-x , 



and integrating, 
^ An 



\ 



(y-y J = wln(xixi } + VnC, 



in(t'(x-x )^] 

/■ 

•Taking antilogs, wc obtain 

y - y, = Cu-x,)^ y * . • 

Let us insert the given condition, which is the 
starting point (x^,y^) I^r the rescue v^essel. Then 



''0 * 



C(x,-x,)^ 
C = 



y . y 

(x -X )^ 



so that finally 



0 - X.J 



Whafkind of path is this? For example, if u = ^, 
this path would be a parabola, leading directly from 
the point {x^.y^) to the point of disaster, (Xj,yj). 
(The vertex ot the parabola is located at (Xj,yj).) 

Model Exam. • - ^ 

^ a) From the form of the surface equation, the lar- * 
^gest value 2 can possibly have occurs when 
X « y « 0;^ indeed 2 = H is the height of the 
mountain! Als6 z - 0 would indicate ground 
level. 'Ihus the family of level curves would 



C «= H - a(x^ ♦ \iy\) , 



0 < C < .H , 



or 



S-l- 



H - C 



« X + uy' 
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This is a family of ellipses, with the ratio of 
the major axis to minor axis remaining constant-. 

b) ^f(x,y) = - 2axi - Zauyl^' 
?f(l,2) = - 2ai - 4awj .' 

c) To find the level\ curve passing through (1,2) » 
insert x = 1, y = 2 into the familv of level 
curves, obtaining (H - C)/a =^ 1 + 4 ; using this 
value for (H - C)/a, the particular member of^ 
the family passing through (1,2) can be identi- 
fied as 



d), 



e) 



1 + 4u 



wy' 



The -tangent vector along a level curve is per- 



pendicular to the Gradient, so from b) , 



dX 



4aVi 
-2aX 



2aj 



dY 
J? 



-2awY 



dY 



dY 

dX 
3? 



-2auy _ wY 

"^ISoc " X 



f) 



dY ^ wdX . , y , ^ ^ ^ ' 

-y " ~T~ ' ulnX + InC ; 



InY InCX^ ; 



Y = CX' 



Inserting the condition, y = Cx^, C = 



thus 









' X(s) 






g) ■ 


'. Y(s): 


Y(s) 


- A 








Z(s) 


= H - 



and 



t> <. s <. X ' 
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1. DESCRIt>TION OF THE PRORI FM 



The orb-it of a satellite aroifnd the earth may 6e 
considered to be determined by the gravitational inter- ■ 
.action between it and the earth alone. Since the sun 
.also affects the motion of the satellite, a small amount 
of error is introduced. However, the earth's gravity is 
the most important gravitational factor 

0 

In the same way, the description of the motion of a 
planet around the sun may be viewed as the result of 
these -two bodies' mutual gravitational attraction. Again 
other planets affect this motion, but the sun is the ma- 
jor influence. on the orbit of a planet. "The dynamics of/ 
the history of the solution of the problem of describing 
• the motion of a planet can be recreated in a short time 
using the modern conveniences of vector differentiation. 
The central ideas and facts are the Inverse Square Law 
and Kepler '.s three Laws. 

The cbmmon thread of the satellite and#anetary 
.motions reappears in modern physics on the sub^molecular 
-level 'in the form of Coulomb" pot^ij'tials . 



2. STATEMENT OF THE LAWS 



« 



We imagine the situation'of two point massesT'one 
•mass very much larger than the other. .The effect of this 
assumption is that we let the position of the larger mass 
be fixed. Put the origin/of the -coordinate system at the 
larger, mass M. For convenience, call it the Sun. The 
motion of the smaller mass m, iiow called a planet de- 
scribes a curve in 3-space. In' Section 6 we will show® " 
that the motions of interest are planar. We win de- ' ^ 
scribe the motions in polar coordinates. The path of the 
pUnet is written in parametric form (r (.t) , e (t) ) -where t 
jis ^time. ^ 

92 - ■ 



FLANET mass ro 



i 




' SU.N mass W 



rnn.Jn!* P^^"'"' """"8 around tlfe sun in a fixed 

orovf^^f K T''"'- 'heliocentric view of the sola^ system ' 

provided by Aristarchus c. 310-230 B.C. and resurrected by 
Copernicus. 1473-1543 led to the discovery of- Kep[er's Laws 
and the Inverse Square. Law. Nepier s Laws . 

- - Two of Newton's Laws are germane to our/discussifn. 
The first is • '. ' ' 



(1) 



mA 



that is, forc'e F gnd acceleration A are v^ctdr quantities 
and they are proportional with the proportionality con- - 
stant being the -mass." In our problem, the acceleratio?? 
that accounts for- the vector motion (r,6) is due* to the 
external forces via equation (1). ' . 



Secondly, we have^ the Inverse Square Law 



(2) 



-GmM 



where G is ^ ^constant which dfepe;ids on, the units' of mea- 
surement, but not on ihe solar Vystem, and u 'is a unit 
vector direct^, from the origin to the mass.M. 

We will discuss three of .Kepler's Laws, . 
Kepler's First Law. .The radius vector to the planet 
sweeps out. area at a constant rate with respect to time, 
that is,the areaof the shaded region is t - t |- 
where .X is a constant, . , ' 
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Figure 2. The radius vector sweeping out area. The 
first of Kepler's A^aws was the .computation of the rate in which 
this area is swept out. It is a constant, even for domets that' 
pass the sun only once. % 



Keplef's gecond Ifaw . The planet's orbit is- ellipse 
j with the sun at a focus. 

Kepler's Third law. If a is the semi-major axis of an 

elliptical orbit and T the time to complete one orbit, ' 
* ■ 

then'T^/a' i$ a solar system constant. In' other words, 
if the quantity T^/a^ is computed for two different plan- 
ets of tfie solar system, it is the s'^ame in both ca^es. 

It^only depends on the mass of the sun and units. - 
f » 

3.* HISTORY OF THE PROBLEM 



ERIC - 



The problem of trying to ex'plain 'tlje p.lanetary mo- 
tion goes back, to antiquity. The part of the history 
which we describe begins in the sixteenth century with 
Copernicifs (1473-1543) who was a proponent of the helio- 
centric view of astronomy, i.e. that the "sun is -at the 
center of things." 

T^cho Brahe (1^546- 1601 > Opposed the Copernicaji theory 
on religious grounds.- It was thought that anything that 



did not put the earth at the center degraded humanity. 
Nevertheless, this ast;*onomer ' s carefuLwork was a major 
contribution in validating the Copernican approach . He 
ha/l received a comi^ission from King Frederich ri of .Den-, 
mark to update astronomical tables. His observatory on 
the island of Hyen-^ODn'tai ned no telescope (it was inven- 
ted in 1609) but he was nevertheless able to record a 
grea^t deal of accurate information. 

This^ accurate information was put to. good use by 
Johann Kepler ( 1 571 - 1 650) , _who w^is Tycho's assistant for 
a short time. Tra^ined as a mathematician he took as his 
task, the study of the orbit ofi^ars. He was an ardervt 
supporter of the CQpern\can theory and his life long am- • 
bition was to find the mystical harmony in the skies. 
His detailed study of Mars led to his publishing his 
first 'two laws in 1609 an^ the third some ten years lateT. 

^ Galileo Galilei ( 1 564- 1642) , ^ivho is well known for 
his^ experiments on particles moving under the influence . 
of gravity,' dismfssed Kepler's astronomy because in intro- 
ducing ellipses he was departing from the more perfect j 
circular motion. Thus the scientist who was to be branded 
a heretic for his scientific vi'^ws rejected, out of hand, 
the work of Kepfer fq^r .very unscientific reasons. 

^n'the.year that Galileo died, another scientist, ' 
Sir Isaac Newton (1642-1727) was born. Newton was well 
acquainted^'wi th the worjc of both Kepler and Galileo and 
of course with ^he iCopernican approach. At the age/of 
25 he discovered that the only g^ravi tational force con- 
sistent witl^ Kepler 's laws ;was the Inverse Square Law, 
He Bid not publish his result -immediately because he at- 
tempted to validate it by doing" calculations* on the orbit 
of the moon. Unfortunately they-^id not ch«ck becaas^ 
some of the data "on distance to the moon was incorrect. 
The cor.rect data i.ndeed did verify the Inverse Square Law. 
He published his result only when it was begun to be pro- ' 
posed by other scientists. '^ O ' 



. An amusing sidelight is that'when asked about the ^ 
posibility o-f the Inverse Square *Law as being correct, 
Newton r^eplied that he had once done the calculations. 
When asked to reproduce nljiem, he could not! Eventually 
he found an error in his second calculation which when 
corrected gave the .correct ansver.* He is generally cre- 
dited witii the discovery* of the Inverse Square 'Law. 

,^In succeeding sections^we try to recreate the scien- 
tific "proces s of going from Tycho's empirical ' data to 
Kepler's Laws tc^Newton's Laws. We will also show tha^ 
starting with Newton's Laws, we can recover Kepler's Laws 

^ 4. MOTIONS described' IN POLAR. COQiiDLNATES 

Suppose that we have a-motion descjribed in polar co- 
ordinates and, R(t) is the position vector. In ordfer. to 
isolate certain aspects of the motion we can introduce^ a' 
local coordinai:e system, as follows. i 



CM 
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• i 
















y R(t) 






• : • « 










^ -.6 =, 0 



Figure 3. The stai^^fard 'unit> vectors for parametric^ 
polar coordinate motions- These form a "moving coordinate 
system" which depend on the position of the particle. This 
modern tool> not available to Kepler and Newton, allows one 
to almost completely dispense with geomepdc and/6r ^rigono- « 
metric arguments. ^ 



*This is a warning to the Scientific neophylie. Keep 'your notebooks 
t)rderly. 



The vectors and U^, are to ^bjp unit^ vectors oriented as * 



indicated. • They 4q' not .depend 'on r bu,t it is easy to see ' 
that U^i^) = \cos 6,sin''e,i and U^O) = (-sin*0,'cos 6), ' 



(3) 



36 ;^Ue 



and 



^ II - 



,Now RCt) = i\(t)Ur(ejt)) is the eq^uation of the motion in 
polar coord ina.t>e? . We-let/v and A be the velocity and 
acce-lerat.iQn vectors, the prime notation means differen- 



- tiation with respect rto time 
V = R' = r'U 



Then 



, ^dUr d6 



r'Up + rU. 



arid 
A 



m 



t ' 

Therefore ^* » 

(4) ' A,= Irw" - r(e')2]u^ + (2r'e'. + re;')U^ 
The coefficients in this vector 

(5) • 



, a^ = r" - r(e')' 



and 
(6) 



a^ = 2r '6', + re'» 




tion respecti 
normal conrponVnt^ 
^motion is on^ a\ circle centered at th^^origin. 



are called t h e^^r^tisi a 1 and angular components of accelei^*^^ 
''These'*are the usual tangential and,''^^^ 
ly in special cases, e.g, when th^*;^ * 



We 'recall al^o that in polar coordinates the aril 
element is rdr so. that the area of the shaded regdcmr^ 
I'n- Figure 2 is given by i 

/7) , S(t)=[d9| rdr.= pir^de - j-^lr^e^^'-,. 

■ - -m^ 



5. DEDUCTION OF THE INVERSE SQUARE LAW 
^. . ^ ^ FROM KEPLER'S LAWS ,> 

Kepler's Laws ar^ emp^irical results based on careful 
observations.^ If one assumes that only forces exterrjal * 
to* thre plane t'' account for its motion, t^ien what mu^t this 
force be? Let us assume Kepler's Laws. 

'* First ob^s^rve tHat ^he fifst law applied to the for- 
mula (7 ) gives , I 
» * * 

• V , • • • • 'I 

^here 'Xlis a constant. 

"The consitant is*divided by 2 to make t'he next eq-ua- 

tion and its further uses simplert Thus 
^ , - • ' ' ' ■ 

' (8) r^e*- =• X . 

If" w% differentiate this relation then . - * ■ 
'^90 > 0 - 2rr'e^, r<^6>' ^ V{2r'^' 'V re*']--' ra' . 

^ ^ WVysee/tJiat* Kepi^'s' firisi£ I^ij.^ = 0, 

,^ th^^t i*s "-J^ie a^ceTei^tion'aivi therefo^-e- .the force is pure- 
^yj^ ly'an-thei-raddai "direction, 'Scujit'^'force^ are cal l^d 



* ' We further assume' that~ the plah6t moyes in a xonic 

sectio-n Ai^eiidix 1), thW-is*^ ■ . • 

"^-..(lO) , ^r(l^+ e cos (6,+ aj) = b\ ^* ) * 

/ Oif fertjnt ia t ipg with respect t-o ,t3^me we get' ' * 

, r' (1- ,e*xos(9 + a)) - re sin (6 ^-^^ a-)^*' v= 0* 

^ I/f we multiply this .equation l?y r and use (8) an^d (10) we 

. ' BrO- Ae 6in(e + a) = 0 . . . '* * • ' 

Now we differentiate again -^,o get ^ |. ' 

. (llj^ / Br" - ,Ae co§(e V"Qi)e' ^ ft- ^, . 
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We solve tlO) f©r,cosX-6 + a) a^d u#e in (11) to get 
Br" - X(| >1)'8' = Q , 4 • ^ 

^ Rearranging, we have 

'\ ^ :i ' 

r B ^ ■ 

On- the left hand side repla.ce X by, (8) and'on the right 

-2 * $ 

sidQ replace 6' Xr by (8) to get: 

, . = r" - r(e')2\= * b7^ 

Since a^ = 0 we have, (compare (1) »and (4)) 

b ' . - 

X^ m V * ' ■ 

p = ^ ^ u ^ the Inverse Square Law, 

Br^ ^ . 

We have shown that t^e - Kepler ' s First and Second' 

Laws imply the Inverse Squar^e Law, , T{ie abov-e calculation 

' ijs no t^ what, Newtdn' $id.' StJmething closer to what he did 

'^is oufjldned, in Exefcise^l where it-5^s shown tjiat th^ Se- 

'cqh^.^nd Third L'aws .'i^ip^'y the ^ Inverse 'Square^ Law , In any 

* . 'ca-se, tihis *Sifetion shows that experimental 'evidence -well 

used c^n lead to nice and {Powerful theoretical results. 
* V • * ♦ 

'- > ■ 1 ■ '-^ —, tt;^ 

Exercise ^ssume^ ^5 Newton did,, that the moon is in a circulat " 

orbit and' thA Keplj^r's Secrond and Third Laws hold. Show that, 0» is 
a constant and therefore it is uniform circufar motion. Introduce 

. the linear speed* v = ds/dt along the circi^e. Show, that a » -(v^/;-). 

. This is the usual formula for oentrij^etal accele^tion. (Jorabine with 
the Thij-d Law to get *thp' Inverse Square Law. 

Exercise 2 . -Justify the- steps in the following withbut doing any 
integrals . " i 

I 2y1>^(l - xVa?) dx '= ab | 2^1 - u^ du = abrr x 

Why. does this show that ^the area o£ ellipse is Trab? We need this 
formula fo'^'fhe area of an ellipse in derivilhg the Third Law. Hint: 
, Evaluate the last integral'by ii^terpretatio^l rather than calculation. 
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KEPLER^S^LAWS AS CONSEOUE NCKS OF NEWTON'S CAWS 



If one decides that F = mA and the Inverse Square 
Law are corxect, do Jepler • s Laws follow? 'This is an im- 
portant question. Kepler's Laws were empirical results'' 
>ased on fanciful hope and data ' wiiich tiad unavoidable? in- 
accuracies. ' Suppose'one can conduct other .e.xperiments 
which verify the Inve^rse Square Law'. If Kepler^s Laws 
are* a consequence o'f the Inverse Square Law^, then Kepler's 
Laws w«ill no-longer be empirical Results from a single 
set o'f datum.A A^physical theory is strengthened by logi- 
cal •impl^catLbns betyeen va'rious empi rical ^sul ts , In 
this section we will >show that Kepler's Laws can be de- 
du(fed from the Inverse -Square Law. , 

. First we give an argument that for any central force 
field, the motions .are pjajiar.^ Suppose thals the planet 

•> is at- a poinl and* that its'moti^ h«s^-a\tangerit vectot' 
T at that point. Draw the plan.e through the sun,, through 
P, and containing the tangent vector T. If Z is Vhft c6*- 
ordiaate normal, to this plane, then the force and hence 
the acceleration in ttie Z direction is zero. That is 

' 2(t J = 0, Z'(XJ - 0 and Z^'(t) = 0 for all t. This im- 
plies Z{t) = 0; the motion is planar. 

r 

Secondly, if the force is a central fx>rce field we 
show Kepler's First Law holds. In the planar polar coor- 
dinates, a^ = 0. Thus we may start at equation (9) and 
retrace our steps to (8) and (7) givi«g the equal ar^a > 
result^. 



Exeroise 3. ^ Do it. 



To get the Second- and Third Laws of Kepler^ we must 
speciiize the central force field to. the 'Inverse* Square 
*L^w. We have » ( ♦ * . \ 



'too- 



(12) . r'' - r(e')2 = - )5M ' * . 

r^ ' ■ • 4r 

and pf course (8). The solution of. these differential 

eq.uations is difficult, but^s made easier by looking ;it 

the g;6aK We want r to be a conic as in, (10). Notice 

^ that, l/r(0> IS very simple. This 'sugg6sts introducing 

the^function w(e) = l/r(e) and deriving a differential 

equation that it satisfie.s. We have 
# * 

- - I dt 1 r' r' ^ 



Then 



d^w 

de^ 



3r[x Jde xe* xe^)[^^^^ ^ "77] «^ 



re^' ^. GM ' GM * 



Then from (8) 



'(13) . d^w ^ _^ GM 



in turns out t^hat every solution of. (13) is given in. the 
form ) . ' 

I 

w(e) = |t ^ ^ cos(e + a) 

for some constants a and B > 0 ' (See. Appendix 2). 
Thus- # 



(14) . . r(8) = 



' B 



e cos(0 + a) ^ 
where ' , , ' 

Hence any such moti*on,'describes a conic section. 
Notice that parabolas a^nd hyperbolas^'are possible. In 
■fact the constants' 3, ;a depend on some initial cotidit^ion. 
If they are selected. so that e > 1, then the mass leaves 
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the solar system. -The only bounded orbits are there- . 
fore ellipses {%nd cities as a special ca^>^This is^ 
Kepler • s* Second Law.^' , " ' • ^ 

^ To derive Kepl^r^s Ihird Law we need tp' compute the 
period X'^nd semi^major^^is a/o£ an eU^ipse. Let b be 
the^-sjemi-jiun^r axis ^arul c;. = " - b^ be/^*he focal length! 

• The-'formula^ for 'r is^^asy. Since^ the area of an 
ellipse is irai), (see Exeixise 2) *we have %y (v) and (8j 



nab 



or 
(15) 




Figure 4*. Aj\ ellipse^with the focus at* the origin. 
This is the^ picture that Kepler saw fh his i?ind*s eye. ' 
.It is one can describe tjeatly in polar coot4inates. * 

The major axis of, an ellipse with ends P and P* has length 

2a*. *The length PF is the mioimum distance from F to the 

ellipse, P*F is, maximum, wjiile the focal distance c from ^ 

the center C is c = 1/2|PP*| - |PF| = a - fPF^ , * ■' . . 
t 

Lo'oking ^t (14) we, see that r is smallest when 
cos(e.+ a) ='1 and largest when cos'(e + a) = -1.. ThUs ^ 



Moreove r 



' (a 



1 + e.^ 



< 



• ft. 

= 5a 



1 +^e 



Using (i6> 



Thus 



B . 



(1 + e)^ 



1 + e 



r 

e) 



^a 



1 + e 

« * 

Putting -this' i nto (15) gives 

■r-l 4tI^ 



[2 - \\ \^)\ 



Ba 



(17) 



GM 



Ba^ - 111 a 

^ut. 4tt^/GM is a Constant that onl^ depends on' the units 

^nd the mas^ of' the sun, a Solar System Constant' Thi,^ 

"IS Kepler»s Third Law. ' ' ^ , ' 
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COMMENTS. 



We have discussed the relationship between Kepler's 
Laws and the ^nvi^rse Square Law in the conj:ext of two 
masses.' We also assumed that the larger' mass 'was fixed. 
In fa'ct, it will wobble- sTightly . What«is fixed is the 
common center o^ mass. When more than two masses are in- 
volved, exact description* oJl the motions can, be ascer.- 
•taineii only in, very" special cases. jrThis* is called the 
n-body problem. It is the focus of a great deal of 
mathematical research.' . 
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8. FURTHER E^RC I S€S ^ 



Exercise 4 . Suppose you observe that" for the earth. a = 1.495 x ioV 
km: and T = 365.25 d^ys. If G = 6.670 x lO"!^ in mVkg-sec^ find ' " 
the mass of the Sun. * ,* ' ' ^ 



Exercise 5 . If the eccentricity of the earth ♦s orbit is given to 
be 0.167322, and a as in Exercise 4 find -the exact equation of fixe 
earth's orl>it. ^ * T . • 



where dG/dr^= -:f. •(You may reoall the formula fo^ ds in polar ccor- 
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Exercise 6 . Yoy know that the usual expression for gravity at sea * 
level is g = -9 . 807 m/sec^ . Use the Inverse^Square, Law and' the fact 
that we .may replace the earth by a pdint mass a't .its center to get ^ 
an exact expression for the Inverse Square. Law with the earth beinc - 
the la;rge mass. Us4 6371 Jem. as the radius of the earth. Find thfe 
mass of the earth. Hint: Compare the twp formulae at the surface 
of the earth. 

Exercise 7 . 15 a' satellite' going abound the earth remains m a cir- 
cular otbit, then cenyjipetal acceleration must balance the acceler- 
ation duetto gravity.* Using this equation," co;npute the linear .speed 
V the satellite must^have if its aljtitude is 100 miles above sea 
level. (Ignore air resistance.) 1 mile = lj50b.35 meters. 

Exercise 8 . The Intelsat seivies geosynchronous satellite jegiains 
^ above a Fixed point on the earth »s equator. ' If it is in aclrrcuiar 
orbit* yoii can deduce its altitude. 'Please do. 

Exercise 9 . OSO 4, launched October 18, 1967 was in an orbit Arith 

radius between 5.37S..an(! 5.697 x lo^m. above s^a level. What is the 

period of its moti|)n? ' * 

* ' ' ' \ 

E^Cercisi 10 , Take the general formula ma^. = f(r)»with 0' removed by 

(8). Multiply by r» and integrate once. TheiT replace X by (8) apd 

derive the formula • - ^ ' 

+ G(r) = constant ' . 



dmates) G is called the potJ^tial energy; j m(ds/dt) ^ is kinetic 
energy. The result of this Exercise is the JConservation* of Ene;rgy. 
Exercise 11 . 'Describe a laboratory experiment that -would show that 
centripetal acceloi-ation'is K(vVr). Newton did this with a thought 
.experiment. Ho Reasoned tJiat; the moon must "fall" in its circular 
orbit in one second the same distance as 'if it ^ere ^dropped from a 
s.tationary position. ' , • ^ ' 

* — • 

^ s 9. ^ HINTS T0:THE SOLUTIONS OF kxERCISES 

Exercise l . From (8) 0» = X/r^ is a constant* For a circle s = r6 

'so ,v = ds/dt; = re' = \/t. Since r" = a a~ = -rfA'l^ = -v?/r. * 
^ ^ « ^ A. 

7rom 2:rr = vT and Keple;-»s Third Law a>. = C/v\ 



Exercise 2. 



The last integral is the area, of a un^ circle. 



Exercise 4 . Solve (17) for M» 

Exercise 5 . Use (16) to compute B. , ■ 

Exercise 6 . Write f(r) = m' C/r^, f (6371. x lO^) - -mg*; 
•Compare with (2) to find ^-y^ 

^Exercige 7 . The equation C/r^ = v^/r can be solved for v in terms 
,of known quantities. 

Exerc'ise 8 » Compute v = ojr, o) the rate of spin of the earth arid ' 
apply th,e equation in Exercise 6 to'find r, or use (17). 

Exercise 9. *- Apply the derivation o^ (16) to fitid a., and thW use 
(17). \ ■ 1 

Exercise 10 . ds^ = dr^ + r^de^ 

Exercise n . Spin a weight at the end of a string. ^ Measure the 

tug of the string. Keeping r fixed, double the speed, remeasure 
the tug.' Keeping v fixed, double r. * . ' 
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\ * ' 10. MODEL £XAM , , 

i^oppose we are considering a p-lanet*s motion around th*e'Sun. • 

a) At what point on trie orb,it is the vector from the Sun to 
the planet turning the most rapidly? f (First Law) 

. r 

b) Recalling that speei v = ds/dt, find 'a formula fer the 

speed of the planet which involves the polar coordinates 

f r and- other constants, but: not 9. (Second and First Laws).^ 
. <4 

On t^ie basis of b) , when is the speed the maximum? ^" 

\, « 

Using^the Third Law> (or otherwise) find the relationship^ be- 
tween'the period of a satellite in a circulBr -orbit *and its 
distance frow the center of the earth/ * ^ 

Describe Why Kepler* s Laws wer^ more acceptable after Newton's 
work.o, ' 

• . * » * ^ 

% . « * / V- 

G^^ven that* the acceleration due.to gravity? at ^ea level on *the 
• earth i» g = -52 ft/sec^. Wnd tthc'vcons^tant in the Inverse 
'square Law in the\ units feet^ and^seconds . Use the'^ radius of 
the. eartiv as 4000 miles apd 5280 feet .'in 1 nvile. . 

Suppose that ^he^ gravitati oiiil ^ fi ^ce wer^' F U^.. 'Woilld it 

- still be true that in a ^tellite* s motion, ""the radius^vector 
sweeps area out i« a constant rate?* . ; . - 
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th6 fo^ci at ^ 



*3 . 



jConrcs m Polar Coord^icnates . If one ebf, 
the ori^gih then the ^qua^ioji of^ any ponic, sectioff ia' pd- 
lar coordinates is of the 'form . 

* ' ^ r - V B • ' 

1 + e cos(e + a) ^ 



7' 

e = 0 
0 < e < 1 
e =*1 

^ "e > 1 



The number e' is called the ecc\entricity and the various 
coiiic sections are given by tjie^ table ^ 

^circle 
e 1 1^ p^s e ' 

parabola ' ' 

hyperbola 

The table is easy 'to remejnber. If e < 1 then the 
denominator is never zero and r is bounde^so we haye an 
ellipse. If » = 1 then the denominator is ?ero precisely 
once as 0 traverses a. complete rotati^^j^ t^fe point '(rvO) 
"jumping across, the open end of'^the parabola at infinity.* 
If e > 1, .then-r is infinite twice "jumping from ©n'e , 
branch -of the hyperbola to the other" each time. 

' . . , APPENDIX 2 • • / * 



Solu tiqns -oi Certain Differentiad Equations . We want to 
ferential < 

" V V 



show 'that the differential equation , ^ * ' 



wj.th v(0) = otg' > v*(0) = has a^ unijc^^r^solut ion . ^ 
Suppose Vj and ^Vj are two solutions and let )\ (#1* for ^ ^ 
function pulled out of a hat) be defined by 

'•h(e) = (V. - v,)2 + (v; - v')2 . ^ • ' 
Then * . ,^ - • ^ ^' * 

^ *^h»(e) = (^^^ - V2)(vj - vj* + ^v; > v»)(vy -'vy) 
. • v;) ijiv^ - v^) ^ (vy - vy)]^ ' 

* • ' • ^ ' . • 16 



= U\ - v»)ll) - D] =^ 0 (Using * twice)"* 

SO h is a constant. But h(0) = 0 means (Vj - v^)^ + 
(v; - = Q jj^yg . Q ^^^^ 

Now one solution of the problem is (just plug it ir\) 



with 



D - a, 

cot a = 



if ^ 0. If a, = 0- take v(e) = D - (D - ajcos 9. v 
Since every solution of the differential equation satis- 
fies the extra conditions at 0 for some and ttj, it is 
of the form v D'+ B c;os(e + a) . To make 6 > 0 replace 
a -by a + TT sinee jcos(e + a + tt) = -cos(e + a) changes 
the sign of ^. • ' * ' 
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DBTAILED SOLUTIONS TQ EXERCISES 

Note that we are using rounded off data so that 'these rijLlts 
.may not agree exactly with published fables of data, • 

Exercise 1. From (8) 9' = X/r^^ is constant. On a circle arclength 
s = er so V i ds/dt = tQ':^ Now r" = 0 so a^. = - r(e»J^-= -(vVr). * 
Compute the circumference of the circle in two ways. Then 2Trr = vT 
and use the* Third Law to write ^ ^ ^ 

r -r ' " r ' jz ^ ^(^1 " "T^ 
Exercise 2 . The first integral is the usual 'area integral for the 
ellipse xVa^ + yVb^ = 1. Let x = au to get the last integral. . 

This integral comEjites the area of the unit circle so is tt.> . 

t " * . ' ^ - 

Exercise 4 . By (17), M = - Put into m, kg, sec. Then 

G = 6.670 X io-»», a = 1,495 x loV m, T = (365'. 25) (24) (3600) sec = 
3.156 X 10' sec. Thus M = 471^ (i .495) ho^ V (6. 670) (3. 156)^ x lo^ = 
1.986 X lo^o kg,. 

Exercise 5. By (16), B = a(l - e^) = (1.495) x io«(l - (.1673)2)^'^ 
= 1.453 X 10^ Now put into (14). ^ * • * 

Exercise 6 . ■f(r) = -(Cm/r^).' At sea .level -(^m/x^) = gm. 
'Writing everyt^ng in m and kg., C* = -gr^ = (9.807)^6.371 x lo^)^ = 
3.981 X io^\ But mC = mMG so M = C/G^ 3.981 x 10^V6.670 x IQ-^^^ 
= 5.968 X lO^^^ kg. 

Exercise .7.. From vVr = C/i^, = C/r, with C as in, Exercise 6. 
But r = (6.371) x lO^ + 100(1609. 35)m, so = [^3.981) x lO^**}/ 
[6.j532 X lo'^], = ,.6095 x lo^ Thus v = 7.807 x lo^ m/sec. 

Exercise 8 . a) Let a> be the- fate of -spin of the earth.. Then 
u) = '2ir/24(3600) =: .727 ^ lO'** rad/dec' Thus v = wr and C/r^ = vVr 
as^in iSxercise 7 IVad to, r^ = C/w^ = 3,9$1 x l0^**/(.727j^' x lO"® = 
75.32 ^^■JO^>(m)\ -Tfien f^'= 4.223' x 10^^ m. ^ * ' 

b) AHemately, using (17) r^ = GMT^/A-n^ and the results of 
Exercise 6 , * • ' 

■■•-\ ■ 18 
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^ = (6.670) X 10*'^* (5.968) x 10^- (24^3600) V4Tr^ 

75.27 ^ 10^^ ^ 
. r*« 4.222 X lo^m. 
Exercise 9 . Look at the derivation of (16).. - 
• a^= ^C^ttax distance + min distance). 

The maximum distance'' is (5.697 x lo^ + 6.371 x .io^)m'and the mmimun 
distance is '(5.375 x lu ' + 6.371 x 10^)m. Thus a =' 6.925*x lo^m. 
Apply (17/. * ' • ' 

T^ = 4TT^aVGM . . 

= 4tt2(6.925 x 10^) V(6.670,x lo" * ^^(3.968 x 10^-) 
= 32.94 X rp^sec^ ^ . 

T= 5.739 X lO^sec. (The listed period is 95.7'rain.) 
Exercise' 10 . From m(r" - r(e')^) = f(r) we get 

so ' . - 

mtr'r". - ^] -«f(r)r« = 0. ' 
Notice that, ? ^ ^ . • 

r 

. and !],''■ ' 

' Ihus ^ ' ' ' ^ . 

This makes the bracket a constant. Now replace X by r^O' to get 

imr(r')^ +. r^CeO^l G(r) = constant''" ^ 

^ • ' . . 

but . ^ • 



so we .get 

' ' - ' . / * 

O J mv^ + 0(r) « constant . . . ^ 



ANSWERS TO MODEL EX^M 



a) From the' first law 6» = A/r^, so 6' is largest when r 
is smallest. 



Now 
r^ 



.1 + e cos I 



so 



dr Be sin 6 



. dS (1 + e cos B)*? 



Thus 



(1 + e cos 6)* 



[(1 + ecos 8)2 +le^ sin^ 6] 



(1 + + 2ecos 6). 



But 



e cos 8 = 



- B 



s6 > 

„2 ^ (e=^- 1)A^ ^ 2Al ^ 
• B2 Br ' 

• I^e^all that B>0, so is largest when r is smallest. 

2. a) The third law is T^ = kr^. dompute the distance travelled 

in one period two way^, vj = 2TTr. Now eliminate T to get 

2 47r^ 1 , . * 

V = — — J or 
k r 

b) Centripetal acceleration mu$t balance the inverse square 
Q acceleration. Thus ^ 



V 

r 



m 

r2 



from which 
„2 ^ GM 



1 ^ 1 
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Newton showed that Kepler's Laws follow from the Inverse Square 
law.. So if one accepts the 'latter then Kepler^s Laws are also 
accg^table. ^ • » > ^ 

If f(r) = -A-^tKen'^ • ^' * 

f (4000- 5280) = -32, . . ' 
thus, , ■ ' 

k = 32(5280)^(4000.)^. ' ^ ' ^ 

Yes, the first law on the. area Hwept out by the radius vector 
holds in any central force field. ' ' ^ 
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